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Abstract 

A methodology is described for the analysis of a transient temperature measure- 
ment made in a flat or curved plate subjected to convective heat transfer, such 
that the surface heat flux, the hot-gas temperature, and the gas heat transfer 
coefficient can be determined. It is shown that if the transient temperature 
measurement is made at a particular point located nearly midway in the thickness 
of the plate there is an important simplification in the data analysis process, in that 
the factor relating the surface heat flux to the measured rate of rise of temperature 
becomes invariant for a Fourier Number above 0.60 and for all values of the 
Biot Number. Parameters are derived, tabulated, and plotted which enable 
straightforward determination of the surface heat flux, the hot-gas temperature, 
and the gas heat transfer coefficient from the measured time, the interior 
temperature of the plate, the rate of rise of temperature, the plate thickness and 
curvature, and the mean thermal properties of tfie plate material at the test 
temperature. It is shown that OFHC copper is the preferred material for 
experimental determination of heat transfer parameters using the described 
methodology. 

The design and method of fabricating and installing a “thermocouple plug” 
which enable measurement of the temperature at an internal position within a 
plate with a minimum of thermal disturbance are described. 

The surface heat flux is determined with an uncertainty the same as the 
uncertainty in the measured rate of rise of temperature and in the mean physical 
properties of the material. The determination of the hot-gas temperature and 
the gas heat transfer coefficient involves an uncertainty which is several times 
the uncertainty of the experimental data. 
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Analysis of Transient Heat Flow to Thick-Walled 

Plates and Cylinders 


I. Introduction 

Many rocket motor and other heat transfer problems 
approximate the conditions of the classical analysis of 
transient heat flow to an infinite slab or cylinder. 

The development of the analytical expression for the 
transient temperature distribution in an infinite slab due 
to the convective heat flux corresponding to a sudden 
change in the temperature of the surroundings is given 
by Kreith (Ref. 1) and other authors. In the nomenclature 
of this report, this expression is 



Jf) 

- 2^ exp (~K F °)^ n 4. sin Sn cos Sn ) 


sin 8„ cos S 


( 1 ) 


assumed to be insulated. The physical properties of the 
plate material ( p , c p , k) are assumed constant over the 
range of temperatures encountered. 

The parameter 8„ is a root of the transcendental 
equation 

cot 8 n = ~ (See Table 1) (2) 

Vi 

INSULATED BACK FACE 

//////////////a y/////////^ (x/L) = 0 


* L 



where a plate of thickness L and uniform temperature 
T 0 is subjected suddenly at t = 0 to convective heat 
transfer on one surface (x/L — 1) from a fluid of tempera- 
ture T g having a boundary heat transfer coefficient h, as 
in Figs. 1 and 2, for a flat plate and a curved plate, 
respectively. The back side of the plate ( x/L — 0) is 


CONSTANT PROPERTIES: p, c , k, h, T 
P 9 

UNIFORM INITIAL TEMPERATURE T Q AT t = 0 

a = Bi =¥- Fo=^- 

P c p k L 2 

Fig. 1. Coordinate system for convective heat transfer 
to flat plate 
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Fig. 2. Coordinate system for convec- 
tive heat transfer to curved plate 


Thus the dimensionless temperature distribution ratio is 
a function of the location within the plate, given by x/L, 
and of the two dimensionless parameters 

_ at 

Fo — jj (3) 


within the heated plate. In the process of developing 
these methods, it was found that if the transient tempera- 
ture measurement is made at a particular location within 
the flat plate ( x/L = 0.57), the heat flux to the surface 
of the plate can be determined independently of the 
value of the Biot Number 



as long as the Fourier Number 



is higher than approximately 0.60. This “direct” deter- 
mination of the heat flux may be all that is required in 
the way of data analysis in many instances, and it also 
simplifies tile procedure for determination of the gas 
driving temperature and heat transfer coefficient. 


and 



where 

k 

. “ P c n 


(4) 


(5) 


Inspection of Eqs, (1-5) shows that, with respect 
to most rocket motor problems, where the initial tem- 
perature and the physical properties are known, the 
temperature— time relationship within the plate is a 
function of two unknowns: the external gas temperature 
T g and the gas heat transfer coefficient h. Therefore, in 
principle, two pieces of information, such as wall tem- 
peratures at two different times or places, or wall temper- 
ature and corresponding rate of change of temperature 
with time at a given time, are sufficient to enable deter- 
mination of those two unknowns and thus permit com- 
plete solution of the problem, as is shown in subsequent 
sections of this report. 


A digital computer solution of transient heat flow to a 
series of thick-walled cylinders enabled determination of 
dimensionless temperature-response factors correspond- 
ing to those obtained from the flat plate analysis. This 
enabled extension of the methodology developed herein 
to the solution of transient heat flow problems involving 
curved plates. 

The charts and procedures presented here can be used 
in the analysis of rocket motor and other heat transfer 
problems which approximate the analytical assumption 
of abrupt initiation of convective heat transfer to a flat 
or curved plate that has constant physical properties and 
that is initially at a known uniform temperature. Separate 
appendices detail ( 1 ) the techniques developed for fabri- 
cation and installation of “thermocouple plugs” enabling 
relatively undisturbed measurement of the temperature 
at a specified point within the interior of a thick plate or 
cylindrical wall, (2) the problems involved in selecting 
metallic materials suitable for use in conducting transient 
heat transfer tests, and (3) the solution of a typical 
transient heat transfer problem. 


In this report, Eq. (I) is used as a basis to develop 
equations, charts, and methods which enable determina- 
tion of the heat flux to the surface of a flat plate, the 
external gas convective heat transfer driving tempera- 
ture, and the gas heat transfer coefficient from transient 
temperature measurements made on the surface of or 


II. Development of Equations for Transient 
Heat Flow in Flat Plates 

The equations used in developing the heat transfer 
analysis procedures described herein are obtained by 
operating on Eq. (1) and by appropriately combining 
the resulting equations with themselves and with the 
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basic equations for convective and conductive heat flux. 
These latter equations are: 

(Convection) q = h(T g — T w ) (6) 

(Conduction) q = +k (-jt ^ (?) 

\ (IX / (x/L) =1.00 

(note that x is negative in proceeding from the heated 
surface to the interior of the plate, using the coordinate 
system shown in Fig. 1 ) . 

A temperature ratio expressed in terms of the tempera- 
ture rise from the initial temperature of the plate can be 
obtained by a transformation of Eq. ( 1) : 


r(T.-r 0 n_ ^ r(^-nm 

_(n-T 0 )J V l(T 0 -T B )if 


« sin 8„ cos 8„ f j- J 

l- 2 2 exp (— 8 ’F°) (S „ + sinl „ C os X) 

71 = ] 


( 8 ) 


Equation (8) can be differentiated with respect to x and 
t respectively to obtain 


It is of interest to see if Eqs. (6-10) can be manip- 
ulated and combined to obtain useful expressions in 
only one independent variable, or expressions which have 
particular properties that make them especially useful 
in the analysis and solution of transient heat transfer 
problems. Two such expressions are presented below. 


First, Eqs. (7), (9), and (10) can be combined to 
produce an expression for a heat flux ratio, which turns 
out to be a function of only Fo, Bi, and the location x/L 
at which the rate of change of temperature with time is 
measured. Thus, 



£s,exp(— 8 S r 0 ) 

n=l _ 

« sin S n cos 8« ( j- j 

X 8 " exp ( ~ S " Fo) (8n + Sin sT^sio 


( 11 ) 


This heat flux ratio eliminates the variable T 9t but only 
at the expense of introducing a new variable, the surface 
heat flux q. This ratio turns out to be useful because of 
a specific property, which will be described later. 



2 (T, 


T 0 ) i-^8«exp(-8’Fo) 


sin S n sin 8» (j-j 
(8„ + sin 8 n cos 8„) 

(9) 



2 (T, 


n-i 


sin 8„ cos 8 n 



(8 n + sin 8»cos 8„) 


( 10 ) 


Equations (8) and (10) are analytical expressions for a 
temperature difference and a temperature derivative 
which can be measured with respect to time. The mea- 
surable quantities are seen to depend functionally on 
(T ff - T 0 ), S„, a, Fo, and x/L. When the definitions of 
these quantities are considered, Eqs. (8) and (10) are 
seen to relate the measurable quantities T s , dTJdt, and 
t to the primary variables T g and h, and to the known 
constants T 0 , p, c p , k, x, and L. 


Finally, as a last step in the formulation of the analyti- 
cal expressions, Eqs. (8) and (10) can be combined to 
obtain the ratio of temperature rise to rate of tempera- 
ture rise. Thus 


2a (Tg - T 0 ) _ 
U / dT x \ 
\dt ) 


2Vexp(-SJFo) 


n= 1 


sin 8„ cos 8„ 

(8„ + sin 8 n cos 8») J 


sin 8* cos 8 r 


j. (j) 

£ eX p { - s;fo) 


( 12 ) 


This ratio does not contain the (T ff — T 0 ) term which 
was present in the original expressions, and it does not 
introduce any new variables. Thus it attains the objective 
of formulating an expression which gives the ratio of two 
measurable quantities in terms of only one independent 
variable, the gas heat transfer coefficient h. 
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III. Discussion of Analytical Results 

The dimensionless functions of Fo, Bi, and x/L given 
by Eqs. (8), (10), (11), and (12) have been calculated 
and are included in Tables 2-27 of this report. 1 The 
dimensionless functions given by Equations (8), (11), 
and (12) have been plotted, based on the above-men- 
tioned tabulated data, and are included in Figs. 3-28 of 
this report. 1 In each case, the tables and graphs cover a 
range of Fourier Number and Biot Number and are 
presented for the heated surface of the plate, for the 
insulated back surface of the plate, and for a location 
near the middle of the plate thickness. 

Figures 9 and 3 give the temperature rise as a fraction 

(T g — T 0 ) at the heated surface and at the insulated 
back surface of a flat plate respectively. The temperature 
ratio is plotted as a function of Fo for several values of 
Bi. Figure 9 shows that the heated surface temperature 
increases with time and that the temperature rises more 
rapidly for plates made of material having a low thermal 
conductivity. 2 Figure 3 shows a similar behavior at the 
back surface of the plate, except that there is a time lag 
before any temperature change occurs at the back surface, 
and at any given time the temperature rise at the back 
surface is much less than the temperature rise at the 
heated surface. Figure 4 shows that the response at a 
position 55% of the thickness of the plate (measured in- 
ward from the back surface of the plate, see Fig. 1), is 
intermediate between the response at the front and back 
surfaces. (The reason for the choice of the 55% position 
will be explained in the next paragraphs.) 

Figures 21 and 14 give the ratio of the heated surface 
heat flux to the rate of change of temperature at the 
heated surface and at the insulated back surface, respec- 
tively, of a flat plate. These curves show that, for a given 
Biot Number, the heat flux ratio does not vary with Fo 
for Fourier Numbers above Fo ~ 0.60. Further, it is 
noted that when referenced to the rate of rise of the 
temperature at the heated surface, the values of the heat 
flux ratio are more than unity, and increase with increas- 
ing Bi, see Fig. 21. 


1 The tables and figures are presented in groups which include 
results from both flat plate and curved plate calculations. The flat 
plate results are identified as corresponding to L/R { = 0. The 
means by which the curved plate data were calculated are de- 
scribed in a later section of this report. 

2 N°te that Fo is a linear function of time, and that Bi is inversely 

proportional to thermal conductivity. 


When the heat flux ratio is referenced to the rate of 
rise of temperature at the back surface of the plate, the 
values of the ratio are less than unity, and decrease with 
increasing Bi, see Fig. 14, This “mirror image” behavior 
led to speculation that there might be a “neutral plane” 
position for temperature measurement such that the cor- 
responding heat flux ratio would be independent of both 
Fo and Bi. Accordingly, the heat flux ratio was calculated 
for several assumed interior reference temperature mea- 
surement positions, and it was found that the heat flux 
ratio for a flat plate did indeed converge on unity for all 
Bi and for Fo > 0.60 when the rate of change of tem- 
perature was measured at a position 57% of the thick- 
ness inward from the back surface ( x/L = 0.57). This is 
shown in Fig. 20. 

For various practical reasons, it was decided to use 
the 55% position ( x/L = 0.55) as a standard location 
for measuring temperatures in the interior of a heated 
wall. 3 Figure 19 shows the heat flux ratio referenced to 
x/L = 0.55; the heat flux ratio is very close to unity at 
low Bi, and has only decreased to 0.98 when Bi = 5.0. 

The dimensionless ratio of the temperature rise' to the 
rate of temperature rise for a flat plate, as given by 
Eq. (11), is plotted on Figs. 28, 23, and 22, corresponding 
to temperature measurements made at the heated surface, 
at x/L = 0.55, and at the insulated back surface, respec- 
tively. These plots are of interest because both of the 
coordinates are dimensionless ratios which can be evalu- 
ated from experimental data and known thermal prop- 
erties. This makes possible the direct determination of 
the Biot Number by interpolation between the lines 
on the graphs once the appropriate temperature response 
data has been measured. And once the value of Bi cor- 
responding to a given set of data has been determined, 
it is an easy task to evaluate the gas heat transfer coeffi- 
cient h and then the effective gas recovery temperature 
T g , and the surface heat flux. 


IV. Transient Heat Flow in Thick-Walled 
Cylinders 

The problem of transient heat flow in a thick-walled 
cylinder is physically similar to that of transient heat flow 


3 The extension to curved plates of the flat plate analysis just de- 
scribed is given in a following section of this report. Since the 
mean location of the “neutral plane” for curved plates having a 
w-ide range of thickness ratios seemed to be near x/L = 0.55, it 
was decided to use this value for the flat plates too, so that the two 
analyses would be on a common basis. 
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in a flat plate, except that the mathematical equations 
are more complicated. Charts giving the temperature rise 
ratio at the surfaces of, and at several positions within, 
a thick curved wall are given in Refs. 1 and 2; however, 
the dimensionless parameters used in the analyses lead- 
ing to these charts are different from those used in the 
JPL work described herein. The analyses and charts are 
presented only for the temperature rise ratios and do not 
extend to the heat flux ratio and the rate of temperature 
rise ratio developed in the JPL work. Thus the material 
in Refs. 1 and 2 is not directly applicable to an attempt 
to extend the methods developed in Sections I— III of this 
report for the analysis of transient heat flow in flat plates 
to the analysis of transient heat flow in thick-walled 
cylinders. 

It was hypothesized, following the flat plate analysis, 
that similar results might be obtained for a thick-walled 
cylinder (representing a section of a combustion cham- 
ber) if the additional curvature parameter L/R » could 
be accounted for in a convenient manner. In order to 
test this hypothesis, and because they had immediate 
need for specific results, Nunz and Dodge at JPL (Ref. 
3) used a thermal network simulation digital computer 
program to compute the temperature distribution and 
heat flux within a series of thick-walled cylinders heated 
from the inside by a hot gas having a constant tempera- 
ture and heat transfer coefficient. The nodal map con- 
sisted of many elements within a pie-shaped segment of 
the cylinder, insulated on both sides. 

The thick-walled cylinder problem was solved for a 
series of curvature parameter values, L/Rj = 0.214, 0.540, 
0.893, and 1.303, and the results were tabulated and 
plotted to show the same parameters and temperature 
response ratios that had been previously developed for 
the fiat plate, which corresponds to the curvature param- 
eter L/Ri = 0. 

It was found that for the thick-walled cylinders there 
was again a unique location ( at the approximate midwall 
position ) where the heat flux was directly related to the 
rate of rise of temperature at the location, and indepen- 
dent of the Biot Number and the Fourier Number, 
except that this time the constant of proportionality was 
a function of the curvature parameter L/R; (Figs. 29 
and 30). 

The location of the "neutral” position within the curved 
wall, where the Biot Number dependence of the heat 
transfer ratio disappeared, varied from x/L = 0.57 at 


L/Ri = 0 to x/L = 0.53 at L/Ri = 1.30. It was arbi- 
trarily decided, as previously mentioned, to use x/L — 
0.55 as a common "neutral position’ for both flat and 
curved plates, thus effectively suppressing what could 
have been a complicating, though minor, extra variable 
in the solution of curved plate transient heat flow prob- 
lems. The heat flux ratio at the x/L = 0.55 neutral posi- 
tion is shown as a function of plate curvature L/Ri, in 
Fig. 30; this figure also shows the minor variation with 
Bi caused by using an arbitrary fixed location for the 
neutral position. That the finite-difference, computer- 
model-derived thick-walled cylinder results are consist- 
ent with the analytically derived flat plate results is 
indicated by the fact that the heat flux ratio varies 
smoothly from the curved plate values to the flat plate 
value on Figs, 30-40. 

The values of the temperature ratio, the heat flux 
ratio, and the temperature rise to rate of temperature 
rise ratio for the curved plates are given in Tables 2-27. 
The corresponding curves are shown in Figs. 3-28. 

In order to assist the user working with transient heat 
flow problems involving curved plates having curvature 
ratios L/R, other than those for which the present calcu- 
lations were performed, a scries of cross plots was made. 
Figures 31-40 give the temperature ratio 

T x To 
t 7 -t 0 

as a function of L/R ; and Bi for fixed values of Fo , with 
the temperature measurement position at x/L — 0.55 and 
1 . 00 . 

The heat flux ratio for temperature measurement at 
x/L = 0.55 is given as a function of L/R { for Fo > 0,60 
and for Bi between 0.05 and 5.0 in Fig. 30. 

The temperature rise to rate of temperature rise ratio 
for temperature measurement at x/L = 0.55 is given vs 
curvature ratio L/Ri for fixed values of Fo in Figs. 41—45. 

V. Solution of Transient Heat Flow Problems 

The analysis and curves developed thus far can be 
applied to practical heat transfer problems which ap- 
proximate the conditions of the analysis. The two types 
of application which are possible are discussed sepa- 
rately below. 
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A. Prediction 

If the hot gas temperature and the gas heat transfer 
coefficient are known, then the temperature of the plate 
and the heat transfer to its surface can be determined 
directly as a function of time. It is presumed that the 
plate is of such a material that over the temperature 
range of interest the thermal properties do not vary much 
from their mean values, that the plate is initially at a 
known uniform temperature, and that the imposition of 
convective heat flux from the hot gas is instantaneous. 

As described previously, Eq. (8) gives the flat plate 
temperature ratio at a given time as a function of T g and 
h, and of the known constants T 0 , p, c p , k , x, and L. Thus, 
with T g and h specified, the temperature ratio at any 
desired time can be calculated. Graphically, the tem- 
perature ratio for the heated surface, x/L = 1.00, of the 
flat plate is read from Fig, 9, and the heated surface tem- 
perature is obtained from this ratio and the other known 
temperatures. The surface heat flux is then calculated 
using Eq. (6). For curved plates, the procedure is exactly 
the same, except that the curvature ratio must be ac- 
counted for by using Figs. 10-13 or Figs. 36-40. 

The prediction procedure described above is useful in 
estimating the transient heat flux and wall temperature 
rise which will be experienced when a flat or curved 
plate of finite thickness is suddenly subjected to convec- 
tive heating. A typical example is the operation of a 
rocket motor with an uncooled combustion chamber, 
primarily for the purpose of determining propellant per- 
formance characteristics. Here the methodology above 
can be used to estimate the allowable test duration {lim- 
ited by the temperature rise of the interior surface of the 
combustion chamber) and to estimate the heat-loss cor- 
rection that should be used in correlating the propellant 
performance measured under the test conditions with 
the ideal thermodynamic performance of the propellant. 

The detailed steps to be followed in making a tem- 
perature and heat flux prediction are summarized for the 
reader’s convenience in Procedure I (Section VI). 

B. Data Analysis 

A useful application of the analysis and charts devel- 
oped herein is in the solution of the inverse problem: the 
determination from measured data (temperature vs time 
within a plate) of the heat transfer rate, the gas heat 
transfer coefficient, and the effective gas recovery tem- 
perature, for the case of instantaneous application of 


convective heat transfer to a plate at an initially known 
uniform temperature. 

Many actual transient heat transfer systems from 
which data may be obtained closely approximate the 
conditions described above; typical is the operation of 
a rocket motor with an uncooled thrust chamber and/or 
nozzle wall. 4 In the actual operation of a rocket motor 
there is usually a short initial period, during which pro- 
pellant flow and combustion chamber pressure are in- 
creasing, before fully developed convective heat transfer 
conditions are established. Thus it is sometimes desirable 
to “refine” the actual test data to establish an effective 
starting time for the fully developed convective heat 
transfer process. After this step, the refined temperature 
vs time data can be used to obtain the corresponding 
values of the heat flux, the gas temperature, and the heat 
transfer coefficient by a procedure which will be devel- 
oped below. 

The analytical solutions for transient heat flow in a flat 
plate, Eqs. (8) and (10), show that the temperature rise 
and the rate of temperature rise at a given time are 
functions of the primary variables T g and h, and of the 
known constants of the system. The inverse problem is 
to solve for T 0 and h as unknowns. By the rules of 
algebra, there must be as many independent equations in 
the unknowns as there are unknowns to be solved for; 
thus two independent equations in T ff and h must be 
written in order to solve for T g and h. These equations 
are obtained, in principle, by substituting two measured 
values of temperature rise, and/or rate of temperature 
rise, at known times and locations into Eq. (8) and/or (10). 

In this particular case, when Eq. (12) is formed as the 
ratio of Eqs. (8) and (10), the variable T g is eliminated; 
thus Eq. (12) is a single equation using both of the pieces 
of experimental data, but involving only h as an un- 
known. The solution of the simultaneous equations, Eqs. 
(8) and (10), or of Eq. (12) alone, would be very difficult, 
so equivalent graphical procedures have been developed, 
based on calculated and plotted values of the tempera- 
ture ratio, the heat flux ratio, and the temperature rise to 
rate of temperature rise ratio. These equivalent graphical 
procedures require as input the same two consistent sets 
of experimental data as would the inverse analytical 
solution. The graphical procedures are equally applicable 
to curved plate problems if the curvature ratio is ac- 
counted for as an additional parameter. 


4 The application of cylindrical curved plate analytical results to a 
conical nozzle is discussed in Appendix C. 
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From the nature of the experimental transient heat 
transfer problem, it is clear that the inverse solution will 
be defined best by data taken at a time sufficiently far 
into the test that the startup perturbation has been "lost’' 
and the temperature rise is high enough to be accurately 
measured. Inspection of the analytical temperature re- 
sponse curves (Figs. 3-13) and the analytical heat flux 
ratio curves (Figs. 14-21) indicates that the data should 
be taken at a time corresponding to Fo > 0.60 and proba- 
bly does not need to be taken at a time later than that 
corresponding to Fo — 2,0, providing that this latter time 
is long compared to the startup perturbation time. (In 
the latter part of Section VII, it is shown that it is desir- 
able to take the data at Fo of 1.5 to 2.0 or more to obtain 
the best accuracy in determining beat transfer coeffi- 
cient and gas temperature.) 

While two temperatures measured at the same place 
at different times or two temperatures at different loca- 
tions measured at the same time would satisfy the need 
for two consistent sets of data, it is clear that a tem- 
perature rise and the corresponding rate of temperature 
rise measured at a single point at a time corresponding 
to 0.60 < Fo < 2.0 would be the optimum data set for 
solution of the inverse problem, both from the standpoint 
of data measurement considerations and because the cri- 
teria requiring a consistent data pair at the latest possible 
point in time is best satisfied. 

For tire reasons above, the discussion of data analysis 
procedures herein will be confined to methods of han- 
dling a data-set composed of temperature rise and rate 
of temperature rise evaluated at the same place and at 
the same instant in time. 

As a practical matter, the temperature at the insulated 
back surface of a heated plate generally has too slow and 
too low a response to provide data acceptable for inverse 
analysis purposes. The temperature at the middle of the 
plate has adequate response and is quite suitable for 
analysis purposes. The temperature at the heated surface 
of the plate would also be satisfactory, except that it is 
difficult to measure accurately, both because the tem- 
perature may rise beyond the limit of calibration of the 
sensing element and, more important, because it is ex- 
tremely difficult to insert a thermal sensor in such a 
manner as to measure a heated-surface temperature with- 
out disturbing the temperature field in the region in 
which the measurement is made. For these reasons, the 
midwall temperature measurement position appears to be 
optimum and is recommended for transient heat flow 


investigations using the methodology and procedures 
described in this report. 5 

In Section III, the simplifications to be realized from 
measuring the temperature at the midplate position 
( x/L = 0.55) were discussed; and as a result of this, 
figures were prepared to enable data analysis when the 
temperature data were obtained at this specific position. 6 

A procedure by which the measured data of tempera- 
ture rise and rate of temperature rise at a given time can 
be used to determine the hot gas temperature and heat 
transfer coefficient, as well as the heat flux to the surface 
of the plate, is described in detail (see Procedure II). 
The charts needed to use this procedure have been 
prepared for the cases where the temperature mea- 
surement is made at the x/L = 0.55 position for flat or 
curved plates, and, additionally, for the cases where the 
temperature is measured at cither surface of a flat plate. 
Values of the thermal properties of the wall material are 
needed to perform these calculations; values for copper, 
nickel, and mild steel are given in Appendix B. 

The surface heat flux is calculated by using the mea- 
sured rate of temperature rise in conjunction with the 
analytically derived value of the heat flux ratio (see 
Eq. (11) and Figs. 14-21 or Fig. 30), where it is to be 
noted that this ratio turns out to be essentially invariant 
with Fo and Bi and to depend only on the curvature 
ratio, when the temperature is measured at the x/L — 0.55 
position within the plate. 

The method of determining the hot-gas temperature 
and gas heat transfer coefficient is based on Eq. (12). It 
was shown previously in the formulation of Eq. (12) that 
the temperature rise to rate of temperature rise ratio is 
independent of T 9 and depends only on the Fourier 
Number and on the Biot Number, which is directly 
related to the gas heat transfer coefficient. Experimental 
values of the temperature rise to rate of temperature rise 
ratio and the Fourier Number can be calculated from the 
known constants and the measured values of the tem- 
perature rise and rate of temperature rise at a particular 


n A method of fabricating and inserting a “thermocouple plug” for 
measuring the transient temperature within a wall is described in 
Appendix A of this report. 

"The methods described in this report can be used to obtain ana- 
lytical solutions and corresponding charts for flat or curved plates 
where the temperature data are obtained at any designated posi- 
tion within the plate. 
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time. Then the Biot Number corresponding to these two 
experimentally measured parameters can be determined 
by referring to a set of curves representing analytical 
solutions of Eq. (12) (see Figs. 22—28 or the equivalent 
Figs. 41—45), 


Once the Biot Number has been determined in this 
manner, the gas heat transfer coefficient is computed 
from Eq. (4), and the gas temperature and the tempera- 
ture of the heated surface of the plate are found by 
reference to plotted solutions of Eq. (8) for the tempera- 
ture ratio (see Figs. 3-13 or the equivalent figure in Figs. 
31-40). Finally, these results can be used to compute a 
surface heat flux, using Eq. (6), which should compare 
closely with the previously obtained heat-flux-ratio value 
of the heat flux. 


The detailed steps to be followed in making the calcu- 
lations above are summarized for the reader’s conveni- 
ence in Procedure II (Section VII). 


VI. Procedure 1: Prediction of Wall 

Temperature and Heat Flux, Given the 
Gas Temperature and Heat Transfer 
Coefficient 

This procedure enables the prediction of wall tempera- 
ture and surface heat flux as a function of time for sud- 
denly initiated convective heat transfer, when the initial 
wall temperature, the thickness and curvature of the 
wall, the physical properties of the wall material, the 
hot gas temperature, and the gas heat transfer coefficient 
are known. The procedure is as follows: 

(1) Calculate Bi = hL/k. 

(2) Calculate corresponding values of time and Fo = 
at/L 2 . When the problem involves a curved plate, 
it will generally be desirable to use the integer 
values of Fo (Fo = 0.60, 0.80, 1.00, 1,20, 1.40), for 
which the L/R t crossplot curves have been pre- 
pared (see Figs. 30-45), and to calculate the corre- 
sponding times. 

(3) Enter Figs. 9-13, or Figs. 36-40, for x/L = 1.00, 
with L/Ri , Fo, and Bi to read the hot- wall tem- 
perature ratio at several values of Fo, where 

r (r. w -r Q )i f (T x - r,) ~\ 


(4) Calculate 


T w = 



( r»-T 0 )i 

r x/L) = 1.00 



(5) Calculate the heat flux: 


q=h (T g - T w ) 


VII. Procedure II: Direct Method of Analyzing 
Transient Wall Temperature Data to 
Determine Gas Temperature, Gas Heat 
Transfer Coefficient, and Surface Heat 
Flux for Convective Heat Transfer 

This procedure can be used to determine gas tempera- 
ture, gas heat transfer coefficient, and surface heat flux 
for convective heat transfer. The procedure requires that 
the temperature be measured within the heated plate at 
a position for which the requisite charts have been pre- 
pared; in this case at x/L = 0.55 for a curved plate and 
at x/L = 0.00, 0.55, or 1.00 for a flat plate. 7 The procedure 
is as follows: 

(1) Prepare the data. 

(a) Plot the measured T x at ( x/L) tes t ~ 0.55 as a 
function of time. 

(b) Plot vs time some parameter which will serve 
as an indication of the state of development of 
the convective heat transfer process. In the 
case of a rocket motor, this could be the cham- 
ber pressure. 

(c) Select an appropriate effective starting time 
for the fully developed convective heat transfer 
process, and a corresponding initial wall tem- 
perature. 

(d) Select a time-space mean wall temperature 
and tabulate the corresponding properties (p, 
c p , and k), which will be assumed constant for 
the data analysis. 

(2) Select a time near the end of the test. Measure T x 
and dT x / dt from a curve faired through the plotted 
data. Compute Fo = at/L 2 , where a — k/pc p and 
t is measured from the effective starting time. The 
selected time must be such that Fo > 0.60. 


7 The calculation steps descrihed here are specifically for tempera- 
ture measurement at the x/L ~ 0.55 position within the wall. The 
special cases of a flat plate with temperature measurement at 
the hot wall, x/L = 1.00, or at the cold wall, x/L = 0, can be 
handled in a manner analogous to that shown. 


8 


JPL TECHNICAL REPORT 32- 1588 



(3) Calculate the heat flux directly from the measured 
rate of temperature rise and the heat flux ratio 
(read from Fig. 30 for the curvature ratio corre- 
sponding to the test conditions when T x is meas- 
ured at x/L — 0.55; Bi is a minor variable on Fig. 
30). Thus 


(6) Enter the temperature ratio curves (Figs. 4-13, or 
the equivalent Figs. 31-40) with L/R,, Fo , and the 
value of Bi determined above and read the tem- 
perature ratio corresponding to x/L = 0.55 and to 
x/L = 1.00 (the heated surface). Then calculate 
T g and T w from the following equations: 



chart, 

(x/L)^0.5 G,(L/Ri) 


(13) 


(4) Compute the dimensionless ratio 


2a ( T x — r 0 ) ' 

u (dr; 
dt 


( 14 ) 


Note: For a curved plate it is best to select the 
time to correspond to the highest integer value of 
Fo for which the analytical curves are crossplotted 
vs the curvature ratio and which falls within the 
test data. 


(5) Enter the temperature rise to rate of temperature 
rise ratio curves for x/L = 0.55 (Figs. 23-27, or the 
equivalent Figs. 41-45), with L/R iy Fo, and the 
value of 


■2a (T x — T„)' 
L 2 ( dTr 
dt 


and interpolate to find the value of Bi which cor- 
responds to these input parameters. 


— . 

(T.- 

•^'o)test, (x/L)= 0.55 

(15) 

>) ~ 

r ( t .- 

ToTl 


L(n- 

Td) Jcliart. (ar/£() = 0.55 


T g = 

= [(T, - 

- T 0 ) + T 0 ] 

(16) 


(T w - To) 


(T x 1 o)te.st 

(*/£-> = 0 


/ 

.bb'| 


[(T x ~ T 0 )1 

) 

_(T,-r 0 )_ 

(«/£.) = 1.00 \ 

~(T. - To)" 


L (Tg - T,)J 

(*/r,) = 0.55 / 


r w = [(T w - To) + To] 


(17) 

(18) 


(7) Compute the heat transfer coefficient from the 
equation 


_ Bi k 

h ~~ (19) 

using the value of Bi obtained in Step 3. 

(8) Compute the surface heat flux from Eq. (20), 
where the quantity in the curved brackets is a 
moderately weak inverse function of Bi. This value 
of heat flux should be reasonably close to the more 
accurate value obtained in Step 3: 



Uncertainties involved in the determination of heat 
flux, gas temperature, and gas heat transfer coefficient 
by the methods of Procedure II are discussed below. 

Heat flux, as computed from Step 3 of Procedure II 
(Eq. 13), has the lowest uncertainty of any of the quanti- 
ties determined. Equation (13) involves only the physical 


properties and dimensions of the plate material, tire 
measured rate of temperature rise, and the analytically 
determined heat flux ratio as read from a prepared chart. 
Each of these quantities can be looked up or measured 
with an uncertainty of a few percent or less. (The heat 
flux ratio is invariant with Fo and is almost invariant 
with Bi.) 
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A potential source of indirect error arises if the wall 
thermocouple is not located at exactly the x/L = 0.55 
position for which the heat flux ratio chart (Fig. 30) is 
prepared. However, examination of Tables 20-22 in the 
text shows that the percentage correction to the mea- 
sured rate of temperature rise for a deviation from the 
desired x/L — 0.55 temperature measurement position is 
of the order of the dimensionless deviation from the 
x/L — 0.55 position and thus should be no more than 1 
or 2% at most, and in a known direction: i.e., if ( x/L) leBl 
< 0.55, the measured rate of temperature rise should be 
increased, and vice versa. 

Heat flux is also computed in Step 8 (Eq. 20) of Proce- 
dure II. Here the evaluation depends on the measured 
temperature rise and on a function of the Biot Number. 
The uncertainty in determining Bi will be discussed 
below, but for purposes of calculating heat flux, it is 
sufficient to note that the function of Bi involved changes 
value only in the order of + 10% for wide variations of 
Bi around a nominal value. Thus a reasonable estimate 
of Bi applied to Eq. (20) will yield a good approximation 
to the true heat flux, though the uncertainty of this value 
will exceed that of the heat flux calculated from Eq. (13). 

The evaluation of heat flux is qualitatively equivalent 
to the evaluation of the product of the gas temperature 
and the gas heat transfer coefficient. It is the separation 
and individual evaluation of these two quantities which 
pose the greatest problem and involve the highest uncer- 
tainties in the data analysis procedure described here. 

Determination of individual values of gas temperature 
and of the gas heat transfer coefficient starts with the 
determination of the Biot Number in Step 5 of Procedure 
II. Inspection of Tables 23-27 and Figs. 2.3-27, as well 
as of the crossplots of Figs. 41-45, shows that the inverse 
determination of Biot Number from the Fourier Number 
and the ratio of temperature rise to rate of temperature 
rise is highly sensitive to small changes in either of these 
latter two quantities. Some calculations of the percentage 
change in Bi corresponding to a percentage change in 
Fo at a fixed value of the ratio of temperature rise to 
rate of temperature rise, and of the percentage change 
in Bi corresponding to a percentage change in the ratio 
of temperature rise to rate of temperature rise at a fixed 
value of Fo, were made for the flat plate case over a 
range of values of Fo and Bi. The Biot Number uncer- 
tainty ratios so obtained are plotted on Figs. 46 and 47. 
Additional calculations, not shown on Figs. 46 and 47, 
indicate that for a curved plate with a curvature ratio 


of L/Ri — 1,3, the Biot Number uncertainty ratios are 
50 to 60% higher than for a flat plate. 

Examination of Figs. 46 and 47 indicates that, for best 
determination of the Biot Number, test conditions should 
be designed so that data are obtained in the region of 
Fo > 1.0 to as high as 2.0 or more if possible and with 
Bi > 0.50. Even in the region of Fo — 2.0 and Bi = 1.0, 
the percentage uncertainty in Bi will be several times 
the uncertainty of the experimental data. 

The uncertainty of the Biot Number determination by 
the methodology described herein is in sharp contrast 
to the uncertainty of the heat flux determination. For 
values of Fourier Number above Fo = 0.60, the heat 
transfer ratio is invariant with Fo and nearly invariant 
with Bi (for the x/L— 0.55 temperature measurement 
position) and depends only on the plate curvature ratio. 
Thus the heat flux uncertainty is only the uncertainty in 
the measured rate of temperature rise and in the average 
values of the physical properties of the material. 

As shown in Step 7 of Procedure II, the gas heat trans- 
fer coefficient is directly related to the Biot Number, and 
any uncertainty in Bi will be directly carried over as an 
uncertainty in the gas heat transfer coefficient. The un- 
certainty in the determination of the hot-gas temperature 
will be of approximately the same magnitude as that in 
the heat transfer coefficient and in the Biot Number. 

Thus, while in principle the methodology developed 
here and summarized in Procedure II enables analysis of 
transient heat flux test data to obtain separate values for 
heat flux, gas temperature, and gas heat transfer coeffi- 
cient, in practice the methodology leads principally to a 
unique means of accurately determining the heat flux, 
while the calculated separate values of gas temperature 
and gas heat transfer coefficient remain subject to an 
appreciable degree of uncertainty. 


VIII. Conclusions 

The methodology and charts presented here enable 
prediction of heat flux and surface temperature for con- 
vective heat transfer to a flat or curved plate, if the hot 
gas temperature and heat transfer coefficient are known. 

The methodology and charts also enable experimental 
determination of the heat flux, the hot-gas temperature, 
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and the gas heat transfer coefficient, providing that the 
temperature at a position midway in the thickness of 
the plate can be measured as a function of time. 

Because of the nature of the analytical solution to the 
convective heat transfer problem, the hot-gas tempera- 
ture and the gas heat transfer coefficient determined 
from the midwall temperature measurement have an 
uncertainty several times that of the experimental data 


involved in their determination. However, the heat flux 
is determined with the same uncertainty as the experi- 
mental data. 

Despite its limitations, the methodology is useful in 
that it provides a new technique for experimentally mea- 
suring convective heat flux and for obtaining a reason- 
able estimate of the corresponding hot gas temperature 
and heat transfer coefficient. 


JPL TECHNICAL REPORT 32-1588 


11 



Definition of Symbols 


Bi Biot Number = hL/h 

c p specific heat of plate material 

Fo Fourier Number = at/L 2 

h gas heat transfer coefficient 

L thickness of flat or curved plate 

q heat flux at heated surface of plate 

Ri radius of curvature of inside surface of cylinder 

T 0 uniform initial temperature of plate 

T g hot-gas recovery temperature 

T w temperature of heated surface of plate 

T x temperature at interior position x within the plate 

x location of position within plate, measured inward 
from the outside (insulated) surface of plate (see 
Fig. 1) 


a thermal diffusivity of plate material = k/pc„ 
8 W root of transcendental equation (see Table 1) 
p density of plate material 


Temperature ratio = 


r (r. - r 0 )l 

_(T g - r 0 )_ 


Heat flux ratio = 



Ratio of temperature rise to rate of temperature rise 


■2a(r,- r 0 )- 
L* / dT z 
dt 
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Table 1. Roots of transcendental equation cot S n = sJBi for flat plate solution 


Bi 


5* 



.05 

.22176 

3.1574 

6.2911 

9.4301 

.10 

.31105 

3.1731 

6.2991 

9.4354 

.15 

.37788 

3.1886 

6.3070 

9.4407 

.20 

.43284 

3.2039 

6.3148 

9.4459 

.25 

.48009 

3.2191 

0.3227 

9.4512 

.30 

.52179 

3.2341 

6.3305 

9.4565 

.35 

.55922 

3.2989 

6.3383 

9.4618 

.40 

.59324 

3.2636 

6.3461 

9.4670 

.45 

.62444 

3.2780 

6.3539 

9.4722 

.50 

.65327 

3.2923 

6.3616 

9.4775 

.55 

.68006 

3.3064 

6.3693 

9.4827 

.60 

.70507 

3.3204 

6.3770 

9.4879 

.65 

.72851 

3.3341 

6.3846 

9.4931 

.70 

.75056 

3.3477 

6.3923 

9.4983 

.75 

.77136 

3.3611 

6.3998 

9.5035 

.80 

.79103 

3.3744 

6.4074 

9.5087 

.85 

.80968 

3.3874 

6.4149 

9.5139 

.90 

.82740 

3.4003 

6.4224 

9.5190 

.95 

,84426 

3.4131 

6.4299 

9.5242 

1.00 

.86033 

3.4256 

6.4373 

9.5293 

2.00 

1.07687 

3.6436 

6.5783 

9.6296 

4.00 

1.26459 

3.9352 

6.8140 

9.8119 

5.00 

1.31384 

4.03354 

6.9096 

9.8928 


h 


*7 


12.5703 

15.7111 

18.8522 

21.9934 

12.5743 

15.7143 

18.8549 

21.9957 

12.5783 

15.7175 

18.8575 

21.9980 

12.5823 

15.7207 

18.8602 

22.0002 

12.5862 

15.7239 

18.8628 

22.0025 

12.5902 

15.7270 

18.8655 

22.0048 

12.5942 

15.7302 

18.8681 

22.0071 

12.5981 

15.7334 

18.8707 

22.0093 

12.6021 

15.7366 

18.8734 

22.0116 

12.6060 

15.7397 

18.8760 

22.0139 

12.6100 

15.7429 

18.8787 

22.0161 

12.6139 

15.7460 

18.8813 

22.0184 

12.6178 

15.7492 

18.8840 

22.0207 

12.6218 

15.7524 

18.8866 

22.0229 

12.6257 

15.7555 

18.8892 

22.0252 

12,6296 

15.7587 

18.8919 

22.0275 

12.6336 

15.7618 

18.8945 

22.0297 

12.6375 

15.7650 

18.8971 

22.0320 

12.6414 

15.7681 

18.8998 

22.0342 

12.6453 

15.7713 

18,9024 

22.0365 

12.7223 

15.8336 

18.9547 

22,0815 

12.8678 

15.9536 

19.0565 

22.1697 

12.9352 

16.0107 

19.1055 

22.2126 


Table 2. Temperature ratio (T* - T 0 )/(T p - T 0 ) vs Fourier and Biot Numbers; L/R { = 0; x/L = 0 


Ratio for indicated Biot Number 


Fo 

Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

.1 




.00227 

.2 





.3 

la 

.01393 

.02702 

.03933 

.6 

.02115 

.04126 

.07862 

.11256 

.8 


.05959 

.11244 

.15951 



.07761 

.14507 

.20404 

1.2 

.04959 

.09529 

.17651 

.24622 

1.4 

.05889 

.11262 

.20680 

.28617 

1.6 


.12963 

.23597 

.32400 

1.8 

.07722 

.14631 

.26407 

.35983 



.16267 

.29113 

.39375 


Bi- 

0.5 

Bi — 
0.7 

Bi = 
1.0 

Bi = 
2.0 

Bi = 
5.0 

.00368 



.01222 

.02269 

.02740 

.03674 

.04936 

.08211 

.13512 

.06185 



.17222 

.26367 

.17175 

.22144 

.28232 

.41241 

.55976 

.23940 

.30429 

.38097 

.53398 

.68827 

.30162 

.37841 

.46614 

.63044 

.77928 

.35875 

.44464 

.53960 

.70694 

.84372 

.41121 

.50381 

.60295 

.76760 

.88934 

.45938 

.55668 

.65758 

.81571 

.92165 

.50361 

.60392 

.70470 

.85386 

.94452 

.54422 

.64612 

.74533 

.88411 

.96072 
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Table 3. Temperature ratio (T x - T 0 )/(T ? - T 0 ) vs Fourier and Biot Numbers; L/R i = 0; x/L = 0.55 


Fo 




Ratio for indicated Biot Number 




Bi - 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi~ 

0.3 

Bi - 
0.5 

Bi = 
0.7 

Bi = 
1.0 

II o 

II o 

.1 



.01417 




.06113 


.17745 

.2 


.01747 

.03392 




.13644 

.21769 

.33329 

.3 

.01387 




.11779 

.15432 

.20081 


.44416 

.6 


.05520 


.14873 

.22448 

.28668 

.36102 

.51239 

.66976 

.8 

■1 

.07331 

.13746 

.19387 

.28795 

.36271 

.44897 

.61335 

.76618 

1.0 

.04732 


.16918 

.23659 

.34621 



.69339 

.83444 

1.2 


.10849 

.19974 



.49129 


.75686 

.88278 

1.4 


.12558 

.22917 

.31536 

.44881 


.64657 


.91700 

1.6 


.14234 

.25752 

.35165 


.59392 



.94123 

1.8 


.15877 

.28482 

.38601 

.53531 

.63719 

.73714 

.87875 

.95839 

2.0 

.09304 

.17490 

.31113 

.41855 

.57332 

.67585 

.77331 

...90385 




Table 4. Temperature ratio ( T x - 

T 0 )/(J g - T 0 ) vs Fourier and Biot Numbers; L/R i 

= 0.214; x/L = 0.55 


Fo 




Ratio for indicated Biot Number 





Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

.6 

.026 

.051 

,095 

.139 

.209 

.268 

.338 

.485 

.643 

.8 

.034 

.067 

.125 

.180 

.267 

.337 

.420 

.582 

.739 


.042 

.084 

.155 

.218 

.322 

.401 

.493 

.662 

.812 

1.2 

.051 

.099 

.192 

.256 

.373 

.461 

.558 

.727 

.862 

1.4 


.115 

.210 

.291 

.420 

.513 

.614 

.779 

.900 

1.6 


.131 

.237 

.326 

.463 

.559 

.662 

.821 

.927 

1.8 

.079 

.148 

.263 

.358 

.502 

.604 

.706 

.857 

.948 


.087 

.161 

.290 

.383 

.539 

.641 

.745 

.886 

.962 



Table 5. 

Temperature ratio (T x - 

T 0 )/<T, - T o> « 

Fourier and Biot Numbers; 

L/Ri = 0.540 

x/L = 0.55 


Fo 




Ratio for indicated Biot Number 




Bi = 

Bi- 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

.6 

.023 

.046 

.086 

.120 

.185 

.239 

.303 

.442 

.598 

.8 


.060 

.112 

.160 

.239 

,304 

.380 

.536 

.695 


.038 

.073 

.137 

.194 

.288 

.362 

.450 

.615 

.768 

1.2 

.045 

.087 

.162 

.228 

,336 

.447 

.510 

.679 

.825 

1.4 

.053 

.100 

.188 

.280 

.379 

.466 

.565 

.734 

,869 

1.6 

.060 

.115 

.210 

.292 

.419 

.512 

.613 

.780 

.901 

1.8 

.069 

.129 

.234 

.321 

.457 

.555 

.655 

.817 

.927 

2.0 

.076 

.142 

.256 

.350 

.494 

.594 

.694 

.851 

.945 
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Table 6. Temperature ratio (T x - T 0 )/(T fl - T 0 ) vs Fourier and Biot Numbers; LfR i = 0.893; x/L - 0.55 


F o 




Ratio for indicated Biot Number 




Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi = 
0.5 

Bi- 

0.7 

Bi = 
1.0 

Bi = 
2.0 

II o 

a* 

.6 

.021 

.030 

.077 

.108 

.169 

.217 

.276 

.408 

.563 

.8 

.028 

.052 

.100 

.141 

.217 

.277 

.348 

.499 

.662 

1.0 

.034 

.066 

.123 

.175 

.261 

.332 

.412 

.575 

.736 

1.2 

.040 

.078 

.146 

.205 

.304 

.383 

.471 

.640 

.794 

1.4 

.047 

.090 

.168 

.235 

.343 

.428 

.522 

.696 

.840 

1.6 

.053 

.101 

.190 

.262 

.382 

.470 

.570 

.742 

.876 

1.8 

(■si 

.114 

.210 

.290 

.417 

.512 

.612 

.782 

.905 

2.0 


.127 

.231 

.316 

.452 

.550 

.652 

.817 

.925 



Table 7. Temperature ratio (J s — 

\)HJ g - T 0 ) VS 

Fourier and Biot Numbers; L/R i 

= 1.303; x/L = 0.55 






Ratio for indicated Biot Number 





Fo 

B i- 

Bi = 

Bi~ 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bl = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

.6 

.019 

.034 

.068 

.098 

.148 

.193 

.250 

.374 

.527 

.8 

.023 

.045 

.089 

.128 

.190 

.248 

.313 

.459 

.6 

1.0 

.030 

.057 

.110 

.157 

.232 

.297 

.372 

.533 

.703 

1.2 

.035 

.068 

.130 

.182 

.272 

.343 

.429 

.598 

.762 

1.4 

.041 

.078 

.150 

,210 

.319 

.387 

.480 

.653 

.810 

1.6 

.047 

.089 

.169 

.237 

.344 

.430 

.526 

.702 

.848 

1.8 

.052 

.100 

.188 

.261 

.379 

.468 

.570 

.742 

.880 

2.0 

.060 

.111 

.206 

.285 

.411 

.504 

.608 

.777 

.907 



Table 8. 

Temperature ratio (7\. 

- T 0 )/(T, - 

T 0 ) vs Fourier and Biot Numbers; L/R t = 

0; x/L - 1.00 


Ratio for indicated Biot Number 

Fo 

Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi - 
0.5 

Bi = 
0.7 

Bi = 
1.0 

Bi = 
2.0 

Bi = 
5.0 

.1 

.01760 

.03470 

.06755 

.09865 

.15610 

.20788 

.27642 

.44639 

.69121 

.2 

.02477 

.04858 

.09354 

.13520 

.20985 

.27459 

.35661 

.54236 

.76847 

.3 


.05938 

.11341 

.16268 

.24899 

.32177 

.41115 

.60160 

.81057 

.0 


.08718 

.16344 

.23046 

.34205 

.43040 

.53173 

.72132 

.88810 

.8 


.10471 

.19427 

.27133 

.39598 

.49120 

.59626 

.77906 

.92078 

1.0 


.12187 

.22391 

.30996 

.44541 

.54542 

.65182 

.82480 

.94391 

1.2 

.07286 

.13870 

.25246 

.34653 

.49078 

.59386 

.69973 



1.4 

.08193 

.15521 

.27995 

.38116 

.53244 

.63714 

.74105 

.88982 

.97188 

1.6 

.09092 

.17140 

.30643 

.41396 

.57069 

.67580 

.77668 

.91263 


1.8 

.09982 

.18728 

.33194 

.44502 

.60582 

.71034 

.80741 

.93072 

,98590 

2.0 

.10863 

.20285 

.35651 

.47443 

.63806 

.74121 

.83391 

.94506 
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Table 9. Temperature ratio (T 2 - T 0 )/< T, - T 0 ) vs Fourier and Biot Numbers; L/R* = 0.214; x/L = 1.00 


Fo 



Ratio for indicated Biot Number 




Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi = 
0.5 

II 

£ 6 

Bi = 
1.0 

Bi = 
2.0 

Bi = 
5.0 

.6 

.042 

.081 

.154 

.214 

.323 

.410 

.510 

.705 

.876 

.8 

.050 

.098 

.181 

.255 

.374 

.466 

.572 

.757 

.910 

1.0 

.059 

.113 

.210 

.291 

.421 

.517 

.627 

.803 

.933 

1.2 

.068 

.130 

.237 

.325 

.467 

.565 

.672 

.842 

.950 

1.4 

.075 

.144 

.262 

.358 

.505 

.608 

.714 

.872 

.964 

1.6 

.083 

.159 

.285 

.389 

.541 

.647 

.750 

.897 

.974 

1.8 

.091 

.174 

.308 

.417 

.575 

.682 

.781 

.927 

.981 

2.0 

.100 

.189 

.332 

.440 

.607 

.711 

.809 

.934 

.988 



Table 10. 

Temperature ratio (T^ - 

T‘ 0 )/<7‘ ir - T 0 ) vs Fourier and Biot Numbers; L/R 

i = 0.540; x/L = 1.00 


Fo 




Ratio for indicated Biot Number 





Bi - 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

.6 

.040 

.075 

.143 

.202 

.303 

.387 

.482 

.673 

.858 

.8 

.047 

.089 

.169 

,236 

.350 

.439 

.540 

.728 

.891 

1.0 

.053 

.103 

.192 

.268 

.392 

.485 

.590 

.774 

.918 

1.2 

.061 

.117 

.214 

.300 

.432 

.529 

.637 

.813 

.938 

1.4 

.068 

.130 

.238 

.329 

.467 

.569 

.677 

.844 

.954 

1.6 

.075 

.143 

.260 

.358 

.503 

.605 

.712 

.870 

.965 

1.8 

.082 

.157 

.282 

.385 

.535 

.640 

.744 

.891 

.973 

2.0 

.090 

.170 

.303 

.412 

.567 

.673 

.773 

.911 

.980 



Table 11. 

Temperature ratio (T ff — 

T 0 )/(T g - T 0 ) vs Fourier and Biot Numbers; L/R 

i - 0.893; x/L = 1.00 


Fo 




Ratio for indicated Biot Number 





Bi = 

Bi = 

fit = 

Bi = 

fit = 

fit = 

Bi - 

fit = 

Bi - 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5,0 

.6 

.037 

.070 

.131 

.187 

.282 

.361 

.455 

.644 

.838 

.8 

.044 

.082 

.154 

.217 

.324 

.410 

.508 

.700 

.874 

1.0 

.050 

.094 

.177 

.246 

.362 

.454 

.556 

.745 

.902 

1.2 

.056 

.106 

.197 

.273 

.400 

.495 

.600 

.783 

.924 

1.4 

.061 

,118 

.218 

.300 

.434 

.532 

.640 

.817 

.941 

1.6 

.069 

.130 

.238 

.327 

.467 

.569 

.674 

.845 

.954 

1.8 

.075 

.142 

.258 

.353 

.498 

.603 

.706 

.869 

.965 

2.0 

.082 

.1.54 

.277 

.377 

.528 

.634 

.737 

.890 

.973 
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Table 12. Temperature ratio (T a - T 0 )J(T g - T 0 ) vs Fourier and Biot Numbers; L/R* = 1.303; x/L = 1.00 


Fo 



Ratio for indicated Biot Number 




Bi ~ 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi = 
0.5 

Bt = 
0.7 

Bi = 
1.0 

Bi = 

2.0 

Bi - 
5.0 

.6 

.032 

.064 

.119 

.172 

.262 

.335 

.427 

.617 

.818 

.8 

.038 

.074 


.199 

.299 

.380 

.477 

.668 

.856 

1.0 

.043 

.085 

.160 

.225 

.334 

.421 

.524 

.714 

.885 

1.2 

.050 


.180 

.250 

.368 

.459 

.566 

.752 

.909 

1.4 

.056 


.199 

.275 

.400 

.494 

.604 

.785 

.928 

1.6 

.061 

.117 

.217 

.299 

.430 

,527 

.639 

.814 

.942 

1.8 

.068 

.128 

.233 

.323 

.460 

.561 

.670 

.839 

.953 

2.0 

.074 

.139 

.250 

.345 

.490 

.592 

.700 

.864 

.964 


Table 13. Heat flux ratio q/\j>c p L(dTJdt)\ vs Fourier and Biot Numbers; L/R t = 0; x/L = 0 


Ratio for indicated Biot Number 


Fo 

Bi = 
0.05 

Bi = 
0.1 

II 

5 6 

Bi = 
0.3 

Bi = 
0.5 

Bi = 
0.7 

Bi = 
1.0 

II c 

CrtS 

° 11 

.1 

3.37053 

3.35259 

3.29027 

3.23913 

3.13227 

3.03505 

2.90450 

2.56361 

2.01977 

.2 

1.13694 

1.35776 

1.33332 

1.31081 

1.26874 

1.23104 

1.18144 

1.06000 

.89044 

.3 

1.10522 

1.09572 

1.07712 

1.05983 

1.02821 

1.00027 

.96415 

.87897 

.76817 

.6 

.99704 

.98902 

.97383 

.95971 

.93430 

.91217 

.88396 

.81915 

.73653 

.8 

.99254 

.98464 

.96970 

.95582 

.93086 

.90913 

.88144 

.81776 

.73616 

1.0 

.99192 

.98405 

.96915 

.95531 

.93044 

.90877 

.88116 

.81764 

.73614 

1.2 

.99184 

.98396 

.96908 

.95525 

.93038 

.90872 

.88113 

.81763 

.73614 

1.4 

.99183 

.98395 

.96907 

.95524 

.93038 

.90872 

.88112 

.81763 

.73614 

1.6 

.99182 

.98395 

.96907 

.95524 

.93038 

.90872 

.88112 

.81763 

.73614 

1,8 

.99182 

.98395 

.96907 

.95524 

.93038 

.90872 

.88112 

.81763 

.73614 

2.0 

.99182 

.98395 

.96907 

.95524 

.93038 

.90872 

.88112 

.81763 

.73614 
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Table 14. Heat flux ratio ql[pc p H 1 6T a l6ty\ vs Fourier and Biot Numbers; L/R i = 0; x/L = 0.55 


Fo 




Ratio for indicated Biot Number 




Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi = 
0.5 

Bi = 
0.7 

Bi = 
1.0 

Bi~ 

2.0 

Bi- 

5.0 

.1 

.92317 




.89864 

.88906 

.87593 

.84125 


.2 

.95746 

.95585 

.95292 


.94516 


.93529 

.92333 

.91364 

.3 


.98208 



.97572 

.97333 

.97052 

.96551 

.96414 

.6 

.99841 

.99769 

.99635 

.99514 

.99302 

.99127 

.98916 

.98494 

.98119 

.8 

.99913 

.99842 

.99708 

.99586 

.99374 

.99196 

.98981 

.98544 

.98140 


.99923 

.99851 

.99718 

.99596 

.99383 

.99204 

.98988 

.98548 

.98141 

1.2 

.99924 

.99853 

.99719 

.99597 

.99384 

.99205 

.98989 

.98548 

.98141 

1.4 

.99925 

.99853 

.99719 

.99597 

.99384 

.99205 

.98989 

.98548 

.98141 

1.6 

.99925 

.99853 

.99719 

.99597 

.99384 

.99205 

.98989 

.98548 


1.8 

.99925 

.99853 

.99719 

.99597 

.99384 

.99205 

.98989 

.98548 

.98141 


.99925 

.99853 

.99719 

.99597 

.99384 

.99205 


.98548 

.98141 



Table 15. 

Heat flux ratio qf[pc p L(dT x /dt)] vs Fourier and Biot Numbers; L/R i 

= 0.214; x/L 

= 0.55 


Fo 




Ratio for indicated Biot Number 





Bi - 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bl = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5,0 

.6 

1.100 

1.100 

1.099 

1.098 

1.097 

1.096 

1.096 

1.095 

1.096 

.8 

1.101 

1.100 

1.100 

1.099 

1.098 

1.097 

1.090 

1.095 

1.096 

1.0 

1.101 

1.100 

1.100 

1.099 

1.098 

1.097 

1.096 

1.095 

1.096 

1.2 

1.101 

1.100 

1.100 

1.099 

1.098 

1.097 

1.096 

1.095 

1.096 

1.4 

1.101 

1.100 

1.100 

1.099 

1.098 

1.097 

1.096 

1.095 

1.096 

1.6 

1.101 

1.100 

1.100 

1,099 

1.098 

1.097 

1.096 

1.095 

1.096 

1.8 

1.101 

1.100 

1.100 

1.099 

1.098 

1.097 

1.096 

1.095 

1.096 

2.0 

1.101 

1.100 

1.100 

1.099 

1.098 

1.097 

1.096 

1.095 

1.096 



Table 16. 

Heat flux ratio q/frc p L(dT x /dt)] vs Fourier and Biot Numbers; LJR i 

= 0.540; x/L 

= 0.55 






Ratio for indicated Biot Number 






Bi ~ 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi~ 

Bi = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

.6 

1.266 

1.266 

1.265 

1.265 

1.265 

1.285 

1.264 

1.265 

1.269 

1.0 

1.267 

1.267 

1.266 

1.266 

1.266 

1.265 

1.265 

1.266 

1.269 

1.2 

1.267 

1.267 

1.266 

1.266 

1.266 

1.265 

1.265 

1.266 

1.269 

1.4 

1.267 

1.267 

1.266 

1.266 

1.266 

1.265 

1.265 

1.266 

1.269 

1.6 

1.267 

1.267 

1.266 

1.266 

1.266 

1.265 

1.265 

1.266 

1.269 

1.8 

1.267 

1.267 

1.266 

1.266 

1.266 

1.265 

1.265 

1.266 

1.269 

2.0 

1.267 

1.267 

1.266 

1.266 

1.266 

1.266 

1.265 

1.266 

1.269 


18 


JPL TECHNICAL REPORT 32-1588 







Table 17. Heat flux ratio q/|>c p L(dT x /dt)3 vs Fourier and Biot Numbers; L/R i = 0.893; x/L = 0.55 


Fo 



Ratio for indicated Biot Number 




Bi = 
0.05 

Bi = 
0.1 

Bi - 
0.2 

Bi = 
0.3 

Bi = 
0.5 

Bi = 
0.7 

Bi = 
1.0 

II o 
« 04 

Bi — 
5.0 

.6 

1.444 

1.444 

1.445 

1.446 

1.447 

1.448 

1.450 

1.455 

1.467 

,8 

1.446 

1.446 

1.447 

1.447 

1.448 

1.449 

1.451 

1.456 

1.467 

1.0 

1.446 

1.446 

1.447 

1.447 

1.448 

1.449 

1.451 

1.456 

1.467 

1.2 

1.446 

1.446 

1.447 

1.447 ' 

1.448 

1.449 

1.451 

1.456 

1.467 

1.4 

1.446 

1.446 

1.447 

1.447 

1.448 

1.449 

1.451 

1.456 

1.467 

1.6 

1,446 

1.446 

1.447 

1.447 

1.448 

1.449 

1.451 

1.456 

1.467 

1.8 

1.446 

1.446 

1.447 

1.447 

1.448 

1.449 

1.451 

1.456 

1.467 

2.0 

1.446 

1.446 

1.447 

1.447 

1.448 

1.449 

1.451 

1.456 

1.467 



Table 18. 

Heat flux ratio q/(>c p L(dT x /dt)] vs Fourier and Biot Numbers; HR i 

= 1.303; x/L 

= 0.55 


Fo 



Ratio for indicated Biot Number 




Bi = 
0.05 

Bi - 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bt = 
0.5 

Bi = 
0.7 

Bi = 
1.0 

Bi = 
2.0 

Bi = 
5.0 

.6 

1.652 

1.652 

1.653 

1.655 

1.656 

1.657 

1.659 

1.666 

1.679 

.8 

1.654 

1.654 

1.655 

1.655 

1.657 

1.658 

1.661 

1.667 

1.679 

1.2 

1.654 

1.654 

1.655 

1,655 

1.657 

1.659 

1.661 

1.667 

1.679 

1.4 

1.654 

1.654 

1.655 

1.655 

1.657 

1.659 

1.661 

1.667 

1.679 

1.6 

1.654 

1.654 

1.655 

1.655 

1.657 

1.659 

1.661 

1.667 

1.679 

1.8 

1.654 

1.654 

1.655 

1.655 

1.657 

1.659 

1.661 

1.667 

1.679 

2.0 

1.654 

1.654 

1.655 

1.655 

1.657 

1.659 

1.661 

1.667 

1.679 



Table 19. Heat flux ratio q/|>c p L(<nydt)] vs Fourier and Biot Numbers; 

HR, = 0; x/L 

= 1.00 






Ratio for indicated Biot Number 




Fo 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 

Bi = 


0.05 

0.1 

0.2 

0.3 

0.5 

0.7 

1.0 

2.0 

5.0 

.1 

.56662 

.57188 

.58374 

.59538 

.61947 

.64413 

.68219 

.81770 

1.28538 

.2 

.79349 

.80435 

.82711 

.85004 

.89728 

.94595 

1.02152 

1.29272 

2.21993 

.3 

.92082 

.93531 

.96504 

.99516 

1.05692 

1.12041 

1.21858 

1.56582 

2.69960 

.6 

1.01135 

1.02812 

1.06202 

1.09632 

1.16602 

1.23703 

1.34559 

1.71999 

2.89390 

.8 

1.01598 

1.03281 

1.06678 

1.10114 

1.17091 

1.24193 

1.35042 

1.72422 

2.89646 

1.0 

1.01662 

1.03345 

1.06742 

1.10177 

1.17152 

1.24252 

1.35096 

1.72460 

2.89660 

1.2 

1.01671 

1.03354 

1.06750 

1.10185 

1.17160 

1.24259 

1.35102 

1.72463 

2.89661 

1.4 

1.01672 

1.03355 

1.06751 

1.10186 

1.17161 

1.24260 

1.35103 

1.72464 

2.89661 

1.6 

1.01672 

1.03355 

1.06751 

1.10186 

1.17161 

1.24260 

1.35103 

1.72464 

2.89661 

1.8 

1.01672 

1.03355 

1.06751 

1.10186 

1.17161 

1.24260 

1.35103 

1.72464 

2.89661 

2.0 

1.01672 

1.03355 

1.06751 

1.10186 

1.17161 

1.24260 

1.35103 

1.72464 

2.89661 


JPL TECHNICAL REPORT 32-1588 


19 







Table 20. Rate of temperature rise (dVdt)/[(2a/L 2 )(T ff - T 0 )] vs Fourier and Biot Numbers; L!R { - 0; x/L = 0 


Ratio for indicated Biot Number 


Fo 

Bi - 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi = 
0.5 

Bi ~ 
0.7 

Bi = 
1.0 

Bi ~ 
2.0 

II o 

a “ 

.1 

.00729 

.01440 

.02834 

.04174 

.06735 

.09135 

.12456 

.21595 

.38221 

.2 

.01780 

.03504 

.06798 

.09896 

.15570 

.20624 

.27229 

.43173 

.65005 

.3 

.02193 

.04292 

.08231 

.11851 

.18260 

.23732 

.30537 

.45325 

.61652 

.6 

.02394 

.04615 

.08590 

.12028 

.17605 

.21856 

.26487 

.34020 

.37982 

.8 

.02382 

.04546 

.08309 

.11435 

.16222 

.19588 

.22902 

.27018 



.02360 

.04462 

.08008 

.10835 

.14901 

.17508 

.19757 

.21428 


1.2 

.02337 

.04376 

.07714 

.10261 

.13683 

.15643 

.17039 

.16992 

.13488 

1.4 

.02314 

.04293 

.07430 

.09718 

.12564 

.13976 

.14694 

.13475 


1.6 

.02291 

.04210 

.07157 

.09202 

.11536 

.12487 

.12672 

.10686 


1.8 

.02269 

.04130 

.06894 

.08715 

.10592 

.11156 

.10929 

.08474 



.02247 

.04051 

.06640 

.08253 

.09725 

.09968 

.09425 

.06720 




Table 21. Rate of temperature rise (dT a ./dt)/[(2a/L 2 )(T ff 

- T 0 )] vs Fourier and Biot Numbers; L/R 4 

= 0; x/L 

= 0.55 

Fo 




Ratio for indicated Biot Number 




Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi ~ 
0.3 

Bi = 
0.5 

II t- 

Bi = 
1.0 

II o 

II o 

S'* 

.1 

IB 

.05245 


.14874 

.23477 

.31184 

.41303 

.65807 

.98215 

.2 

HEggfl 

.04977 

.09512 

.13653 

.20900 

.26987 

.34395 

.49564 

.63354 

.3 




.12834 

.19242 

.24389 

.30337 

.41262 

.49120 

.6 

.02391 

.04575 

.08396 

.11600 

.16564 

.20112 

.23670 

.28294 

.28512 

.8 

.02366 

.04483 

.08081 

.10976 

.15196 

.17953 

.20395 

.22420 

.20181 

1.0 

.02343 

.04397 

.07783 

.10393 

.13951 

.16038 

.17587 

.17778 

.14289 

1.2 

.02320 

.04313 

.07496 

.09842 

.12809 

.14329 

.15167 

.14098 

.10117 

1.4 

.02297 


.07221 

.09320 

.11761 

.12802 

.13080 

.11180 

.07164 

1.6 


.04149 

.06955 

.08826 

.10799 

.11438 

.11280 

.08866 

.05072 

1.8 

.02252 

.04070 

.06699 

.08358 

.09916 

.10219 

.09728 

.07034 

.03591 

2.0 

.02230 

.03992 

.06453 

.07915 

.09104 

.09130 

.08389 

.05575 

.02543 
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Table 22. Rate of temperature rise (dTJ dt/ [2a/ L^)(T g - T 0 )] vs Fourier and Biot Numbers; L/R i - 0; x/L - 1.00 


Fo 




Ratio for indicated Biot Number 




Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

BI- 

OS 

Bi = 
0.5 

Bi = 
0.7 

Bi = 
1.0 

Bi - 
2.0 

II o 

.1 

.04334 

.08439 

.15974 

.22709 

.34057 

.43041 

.53033 

.67703 

.60058 

.2 

.03073 

.05914 



.22015 


.31492 



.3 



.09187 

.12621 

.17764 

.21187 

.24161 

.25443 

.17543 

.6 

.02360 

■s 

.07877 


.14107 

.16116 

.17400 

.16202 


.8 

.02327 

.04334 

.07553 

.09926 

.12896 

.14339 

.14949 

.12814 


1.0 

.02302 


.07271 

.09394 

.11835 

.12805 

.12886 



1.2 


.04167 


.08896 

.10866 


.11113 


.03428 

1.4 

.02257 

■SB 

.06745 

.08424 

.09977 


.09583 

.06388 

.02427 

1.6 

.02235 


.06497 

.07978 

.09161 




.01719 

1.8 

.02213 


.06258 

.07555 

.08411 


KB 


.01217 

2.0 

.02192 

111 

.06028 

.07155 

.07723 


.06147 

.03186 

.00862 


Table 23. Ratio of temperature rise to rate of temperature rise (2o/L z )[(T ir - T 0 )/(dT,/dt>3 vs Fourier and Biot Numbers; 

L/R i = 0; x/L = 0.55 


Ratio for indicated Biot Number 


Fo ' 

Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

.1 

.1372 

.1377 

.1390 

.1401 

.2 

.3484 

.3511 

.3565 

.3618 

.3 

.5627 

.5684 

.5796 

.5907 

.6 

1.1874 

1.2067 

1.2448 

1.2822 

.8 

1.6021 

1.6352 

1.7010 

1.7664 

1.0 

2.0201 

2.0712 

2.1738 

2.2765 

1.2 

2.4421 

2.5156 

2.6644 

2.8150 

1.4 

2.8682 

2.9867 

3.1737 

3.3837 

1.6 

3.2985 

3.4306 

3.7025 

3.9841 

1.8 

3.7331 

3.9015 

4.2514 

4.6182 

2.0 

4.1720 

4.3816 

4.8214 

5.2877 


li 1C 
w ° 

II t> 

a d 

Bi ~ 
1.0 

Bi - 
2.0 

II o 

a* 

.1424 

.1447 

.1480 

.1580 

.1807 

.3722 

.3823 

.3967 

.4392 

.5261 

.6121 

.6327 

.6620 

.7457 

,9042 

1.3552 

1.42.54 

1.5252 

1.8110 

2.3491 

1.8950 

2.0204 

2.2014 

2.7357 

3.7965 

2.4816 

2.6850 

2.9840 

3.9002 

5.8397 

3.1204 

3,4287 

3.8913 

5.3685 

8.7253 

3.8160 

4.2610 

4.9433 

7.2200 

12.801 

4.5736 

5.1926 

6.1632 

9.5549 

18.556 

5.3986 

6.2352 

7.5776 

12.499 

26.685 

6.2972 

7.4022 

9.2178 

16.212 

38.166 


Table 24. Ratio of temperature rise to rate of temperature rise (2a/L 2 )[(Tj. — T (J )/(dT x /dt )] vs Fourier and Biot Numbers; 

L/Ri = 0.214; x/L = 0.55 


Fo 



Ratio for indicated Biot Number 




Bi = 
0.05 

Bi- 

0.1 

Bi = 
0.2 

Bi- 

0.3 

Bi = 
0.5 

Bi = 
0.7 

II o 

a ^ 

Bi = 
2.0 

Bi = 
5.0 

.6 

1.21 

1.24 

1.26 

1.28 

1.36 

1.42 

1.52 

1.78 

2.27 

.8 

1.61 

1.64 

1.71 

1.77 

1.88 

2.00 

2.16 

2.64 

3.57 

1.0 

2.02 

2.08 

2.17 

2.26 

2.46 

2.63 

2.92 

3.72 

5.38 

1.2 

2.43 

2.51 

2.66 

2.80 

3.17 

3.35 

3.76 

5.04 

7.93 

1.4 

2.87 

2.96 

3.15 

3.36 

3.72 

4.12 

4.72 

6.67 

11.30 

1.6 

3.29 

3.43 

3.68 

3.94 

4.44 

4.98 

5.81 

8.67 


1.8 

3.72 

3.88 

4.21 

4.55 

5.23 

5.93 

7.08 

11.17 


2.0 

4.17 

4.37 

4.76 

5.18 

6.07 

7.00 

8.56 

14.27 
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Table 25. Ratio of temperature rise to rate of temperature rise <2a/L 2 )[(T x - T 0 )/(«/T J! /<ft)] vs Fourier and Biot Numbers; 

L/Ri = 0.540; x/L = 0.55 


Ratio for indicated Biot Number 


Fo 

Bi = 
0.05 

Bi = 
0.1 

Bi - 
0.2 

II CO 

3 d 

Bi = 
0.5 

Bi = 
0.7 

II o 

3 ^ 

II o 

3 04 

II o 

S" 5 

.6 

1.21 

1.22 

1.26 

1.28 

1.34 

1.41 

1.48 

1.70 

2.12 

.8 

1.62 

1.64 

1.70 

1.75 

1.86 

1.95 

2.09 

2.50 

3.27 

1.0 

2.02 

2.07 

2.15 

2.24 

2.40 

2.56 

2.78 

3.44 

4.82 

1.2 

2.45 

2.52 

2.62 

2.75 

2.98 

3.19 

3.55 

4.58 

6.83 

1.4 

2.87 

2.96 

3.13 

3.28 

3.60 

3.93 

4.40 

5.97 

9.54 

1.6 

3.29 

3.41 

3.62 

3.84 

4.26 

4.73 

5.38 

7.66 

13.04 

1.8 

3.72 

3.86 

4.14 

4.42 

4.98 

5.60 

6.52 

9.64 


2.0 

4.15 

4.31 

4.67 

5.02 

5.74 

6.52 

7.74 

12.07 



Table 26. 

Ratio of temperature rise to rate of temperature rise (2a/L 2 )[(T 3 — 

L/R i = 0.893; x/L = 0.55 

T 0 )/(dT x /dt>] vs Fourier and Blot Numbers; 

Ratio for indicated Biot Number 

ro 

Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi = 
0.3 

Bi = 
0.5 

II r- 
3 ® 

Bi = 
1.0 

II o 

3 

Bi = 
5.0 

.6 

1.22 

1.24 

1.26 

1.29 

1.35 

1.40 

1.48 

1.68 

2.04 

,8 

1.63 

1,66 

1.70 

1.75 

1.84 

1.93 

2.06 

2.42 

3.10 

1.0 

2.05 

2.09 

2.15 

2.23 

2.38 

2.50 

2.70 

3.29 

4.45 

1.2 

2.47 

2.52 

2.61 

2.71 

2.94 

3.12 

3.43 

4.32 

6.19 

1.4 

2.90 

2.97 

3.10 

3.22 

3.52 

3.79 

4.23 

5.54 

8.38 

1.6 

3.31 

3.42 

3.58 

3.76 

4.15 

4.54 

5.12 

6.96 

11.32 

1.8 

3.74 

3.87 

4.10 

4.32 

4.81 

5.31 

6.10 

8.65 


2.0 

4.15 

4.32 

4.60 

4.90 

5.53 

6.16 

7.19 

10.66 



Table 27. 

Ratio of temperature rise to rate of temperature rise (2a/L 2 )[(7 a , - 

L/R; = 1.303; x/L = 0.55 

T 0 )/(dT x /dt)] vs Fourier and Biot Numbers; 

Ratio for indicated Biot Number 


Bi = 
0.05 

Bi = 
0.1 

Bi = 
0.2 

Bi — 
0.3 

Bi = 
0.5 

Bi = 
0.7 

Bi - 
1.0 

II o 
85 « 

Bi = 
5.0 

.6 

1.22 

1.24 

1.27 

1.29 

1.31 

1.38 

1.45 

1.62 

1.96 

.8 

1.63 

1.66 

1.70 

1.74 

1.77 

1.90 

2.02 

2.32 

2.90 

1.0 

2.05 

2.10 

2.13 

2.20 

2.32 

2.45 

2.62 

3.10 

4.09 

1.2 

2.45 

2.50 

2.60 

2.68 

2.87 

3.04 

3.30 

4.05 

5.58 

1.4 

2.87 

2.92 

3.05 

3.18 

3.43 

3.67 

4.03 

5.12 

4.74 

1.6 

3.28 

3.37 

3.54 

3.70 

4.03 

4.35 

4.84 

6.36 

9.86 

1.8 

3.72 

3.82 

4.04 

4.23 

4.65 

5.08 

5.73 

7.82 

12.87 

2.0 

4.16 

4.28 

4.54 

4.78 

5.34 

5.87 

6.73 

9.50 
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TEMPERATURE RATIO (T - T n )//T - T„) TEMPERATURE RATIO 



Fig. 3. Temperature ratio vs Fourier and Biot Numbers; LfR i = 0; x/L = 0 



Fig. 4. Temperature ratio vs Fourier and Biot Numbers; L/R { =0; x/L = 0.55 
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TEMPERATURE RATIO 


x/L - 0.55 


BIOT NUMBER 



FOURIER NUMBER 

Fig. 5. Temperature ratio vs Fourier and Biot Numbers; L/R 4 = 0.214; x/L = 0.55 



Fig. 6. Temperature ratio vs Fourier and Biot Numbers; L/R { = 0.540; x/L = 0.55 
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Fig. 7. Temperature ratio vs Fourier and Biot Numbers; L/R i = 0.893; x/L = 0.55 
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Fig. 8. Temperature ratio vs Fourier and Biot Numbers; L/R i = 1.303; x/L = 0.55 
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Fig. 11. Temperature ratio vs Fourier and Biot Numbers; L/R i = 0.540; x/L ~ 1.00 



Fig. 12. Temperature ratio vs Fourier and Biot Numbers; L/R i = 0.893; x/L = 1.00 
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HEAT FLUX RATIO TEMPERATURE RATIO 



Fig. 13. Temperature ratio vs Fourier and Biot Numbers; L/R i = 1.303; x/L = 1.00 



Fig. 14. Heat flux ratio vs Fourier and Biot Numbers; L/R i = 0; x/L = 0 
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Fig. 15. Heat flux ratio vs Fourier and Blot Numbers; L/R { = 0; x/L = 0.55 



Fig. 16. Heat flux ratio vs Fourier and Biot Numbers; L/R i = 0.214; x/L = 0.55 
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Fig. 19. Heat flux ratio vs Fourier and Biot Numbers; L/Rj = 1.303; x/L = 0.55 
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Fig. 20. Heat flux ratio vs Fourier and Biot Numbers; L/R i = 0; x/L = 0.57 
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HEAT FLUX RATIO 



Fig. 21. Heat flux ratio vs Fourier and Biot Numbers; L/R i = 0; x/L = 1.00 
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Fig. 22. Ratio of temperature rise to rate of temperature rise vs 
Fourier and Biot Numbers; L/R i = 0; x/L = 0 


Fig. 23. Ratio of temperature rise to rate of temperature 
rise vs Fourier and Biot Numbers; L/R i = 0; x/L = 0.55 
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RATIO OF TEMPERATURE RISE TO RATE OF TEMPERATURE RISE 









RATIO OF TEMPERATURE RISE TO RATE OF TEMPERATURE RISE 



Fig. 25. Ratio of temperature rise to rate of temperature rise vs Fourier and Biot Numbers; L/R i = 0.540; x/L - 0.55 
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RATIO OF TEMPERATURE RISE TO RATE OF TEMPERATURE RISE 



Fig. 26. Ratio of temperature rise to rate of temperature rise vs Fourier and Biot Numbers; L/R { = 0.893; x/L = 0.55 
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RATIO OF TEMPERATURE RISE TO RATE OF TEMPERATURE RISE 
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Fig. 27. Ratio of temperature rise to rate of temperature rise vs Fourier and Biot Numbers; L/R i = 1.303; x/L = 0.55 
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RATIO OF TEMPERATURE RISE TO RATE OF 



Fig. 28. Ratio of temperature rise to rate of temperature rise vs 
Fourier and Biot Numbers; L/R i = 0 ; x/L = 1.00 



Fig. 29. Variation of heat flux ratio with Biot Number and 
x/L value for L/R i = 1.303 
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ig. 30. Variation of heat flux ratio with LfR i and Biot Number 
for x/L = 0.55 and for a Fourier Number > 0.60 


Fig. 32. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 0.80; x/L = 0.55 



Fig. 31. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 0.60; x/L = 0.55 


Fig. 33. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 1.00; x/L = 0.55 
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Fig. 34. Crosspiot of temperature ratio vs curvature ratio and 
Biot Number Fourier Number = 1.20; x/L = 0.55 


Fig. 36. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 0.60; x/L = 1.00 



Fig. 35. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 1.40; x/L = 0.55 



Fig. 37. Crossplot if temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 0.80; x/L = 1.00 
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38. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 1.00; x/L = 1.00 
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Fig. 39. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 1.20; x/L = 1.00 
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Fig. 40. Crossplot of temperature ratio vs curvature ratio and 
Biot Number; Fourier Number = 1.40; x/L = 1.00 



CURVATURE RATIO y'R. 


Fig. 41. Crossplot of ratio of temperature rise to rate of tempera- 
ture rise vs curvature ratio and Biot Number; Fourier Number 
= 0.60; x/L = 0.55 
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Fig. 42. Crossplot of ratio of temperature rise to rate of tempera- 
ture rise vs curvature ratio and Biot Number; Fourier Number 
= 0.80; x/L = 0.55 



Fig. 44. Crossplot of ratio of temperature rise to rate of tempera- 
ture rise vs curvature ratio and Biot Number; Fourier Number 
= 1.20; x/L = 0.55 
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Fig. 43. Crossplot of ratio of temperature rise to rate of tempera- 
ture rise vs curvature ratio and Biot Number Fourier Number 
= 1.00; x/L = 0.55 
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Fig. 45. Crosspiot of ratio of temperature rise to rate of tempera- 
ture rise vs curvature ratio and Biot Number Fourier Number 
= 1.40; x/L = 0.55 
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Fig. 46. Biot Number uncertainty ratio for uncertainty in 
Fourier Number only 
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Fig. 47. Biot Number uncertainty ratio for uncertainty in ratio of 
temperature rise to rate of temperature rise only 
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Appendix A 

Thermocouple Plug Design, Fabrication, and Installation 


I. Introduction 

This report describes the development and application 
of a methodology for analyzing transient temperature 
data obtained during convective heating of a flat or 
curved plate, so as to determine the temperature and the 
heat transfer coefficient of the hot gas. The methodology 
is in its simplest and most effective form when the tem- 
perature measurement is made near the center of the 
thickness of the plate. In order to make use of the 
methodology, then, it was necessary to develop a means 
of measuring the temperature within a metal plate, with 
a negligible local perturbation of the temperature field. 

The technique developed was to fabricate a “thermo- 
couple plug” of the same material as the plate and then 
to braze the plug into a precisely machined and closely 
fitting hole in the plate. Separate fabrication of the 
thermocouple plug enabled special care to be taken in 
installing the fine thermocouple wires and in forming the 
thermocouple junction, using a laser welding process, and 
also enabled inspection of each completed thermocouple 
plug before it was brazed into the plate. 

Brazing of the completed thermocouple plugs into 
machined holes in the plate, with a thin but complete 
metallic bond between the plug and the plate, minimized 
the thermal disturbance due to the introduction of the 
plug into the plates. Additionally, the depth of the finite- 
sized machined hole in the plate could be accurately 
measured before brazing to establish the exact position 
at which the thermocouple junction would be located 
after the brazing was completed. 

The critical steps in the fabrication and installation 
of the thermocouple plugs are described below. 

II. Thermocouple Plugs 

The thermocouple plug was made of the same material 
as the plate; all of the present work was done with cop- 
per. 8 An outside diameter of 0.6-35 cm (1/4 in.) was chosen 
for ease in fabrication and handling. The inside of the 


8 Tellurinm copper was used for the plugs to improve machinability, 

while OFHC copper was used for the plate. The difference in 
thermal properties is assumed to be negligible. 


plug was bored with holes of gradually decreasing diam- 
eter to leave a thin closure at the thermocouple end of 
the plugs as shown on Fig. A-l. The critical part of the 
thermocouple plug was at the closed end, where a dia- 
phragm only 0.0127 cm (0.005 in.) thick remained at the 
bottom of the final 0.09 14- cm (0.036-in.) diameter bored 
hole. The final step in machining the body of the plug 
was to drill two 0.0153-cm (0.006-in.) diameter holes, 
spaced 0.254 cm (0.010 in.) apart, located symmetrically 
about the center line of this thin diaphragm. 

The preparation of the thermocouple wires used in 
the plug was also very critical. (Most thermocouples 
were made with Chromel-Alumel wires, but higher- 
temperature-capability Platinel-Platinel thermocouple 
wire was used for some samples.) Thermocouple wires 
0.0127 cm (0.005 in.) in diameter (.36 gage) were used, 
The first step was to cut the wires into 15.2-cm (6-in.) 
lengths and to form a perfectly spherical ball 0.02033 to 
0.0254 cm (0.008 to 0.010 in.) in diameter on one end of 
each wire; this was done with a controlled pulse of a 
laser beam that was focused on the end of die wire (the 
axis of the wire was aligned with the axis of the laser 
beam). 9 The end of the wire was protected from oxida- 
tion with an argon purge. 

The next step was to thread a pair of thermocouple 
wires through the holes in the diaphragm at the closed 
end of the thermocouple plug and through the similarly 
spaced holes in a piece of 0.0794-cm ( 1/32-in.) diameter 
double-bore ceramic insulating tube, as shown in Fig. 
A-2. The thermocouple plug was then mounted vertically, 
with its closed end up, under the laser. The ceramic tube 
was restrained so that its upper end was about 0.158 cm 
(1/16 in.) below the bottom surface of the diaphragm, and 
lightweight clamps were attached to the bottom of each 
wire to provide a slight tension to keep the spherical 
ball at the top end of each wire in firm contact with the 
top surface of the diaphragm. Then the laser beam was 
focused to cover both spherical wire ends, and a con- 
trolled pulse fused both thermocouple wires into the 
diaphragm to form an integral thermocouple junction. 
The fused end was then flattened by lightly stroking it 
on a sheet of No. 400 wet-or-dry abrasive grit, with the 
plug suitably supported to keep the end square. 


9 The laser used was a CORAJD 150-mA machine. 
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Fig. A-l. Machined thermocouple plug body (tellurium copper) 



Fig. A-2. Assembly of parts for fusing thermocouple junction 
in thermocouple plug 


After the completion of the laser welding and cleanup 
operation described above, the thermocouple plug was 
turned open-end-up, and a high-temperature ceramic in- 
sulating cement was used to cement in the double-bore 
ceramic tubing and part of the space above it; care was 
taken in this last step to keep the thermocouple wire 


leads adequately separated. 10 When the cement had dried 
and set, it was cured at the recommended temperature. 
The remainder of the cavity in the top of the thermo- 
couple plug was filled with a lower-temperature cement 
after the plug had been brazed into the plate. At this 
time, too, woven fiber glass sleeving was slipped over the 
thermocouple lead wires, to be bonded in place by the 
cement. 11 The final configuration of the assembled plug 
was shown in Fig. A-3. 

III. Assembly by Brazing 

A considerable number of copper test specimens were 
used in developing a satisfactory brazing procedure. The 
objective of the brazing operation was to bond the ther- 
mocouple plug to the walls of the hole machined into 
the plate so that there would be an absolute minimum 
of thermal disturbance; i.e., the bond would be con- 
tinuous and thin and would be formed from a high- 
conductivity brazing alloy. These objectives were 
achieved by adhering to the following step-by-step 
procedure: 

(1) The hole in the plate was machined to receive the 
thermocouple plug, and a brazing feed-tube hole 


10 The cement used was Ecoceram 21, made by Emmerson and 
Cumming, which will withstand temperatures of np to 1368 K 
(2000°F). 

n The cement was Englehard Industries CA9R cement, for 812 K 
(1000°F) service. 
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Fig. A-3. Completed thermocouple plug 


was also machined, as shown on Fig. A-4. The 
feed-tube hole was positioned so that there was 
gravity feed of the brazing alloy when the assembly 
was loaded into the brazing furnace. The plug-hole 
was precisely machined with a flat, square bottom 
and a diameter only 0.00127 to 0.00508 cm (0.0005 
to 0.0020 in,) more than that of the plug, depending 
on the brazing alloy used. 

The depth of the hole was such that after assembly 
the effective thermocouple junction was at the de- 
sired position (usually 0.55 of the plate thickness in 
from the unheated surface). Somewhat arbitrarily, 
the effective thermocouple junction was assumed 
to be at the midthickness of the copper diaphragm 
at the tip of the thermocouple plug, and it was 
assumed that there was a 0.00254 cm (0,001 in.) 
thick layer of brazing alloy between the bottom of 
the plug and the bottom of the hole in the plate; 
thus, the effective thermocouple junction was 
0.00889 cm (0.0035 in.) from the bottom of the 
machined hole. 

(2) The bottom of the machined thermocouple plug 
hole was coined with a square-ended punch made 
of hardened steel, This assured that the thermo- 
couple plug could seat firmly on the bottom of the 
hole and would not ride on the corner radius or any 
irregularity on the bottom surface of the hole. 

(3) The thermocouple plug, the hole in the plate, the 
copper feed-tube, and the brazing alloy were thor- 
oughly cleaned. (The cleaning process found to be 
satisfactory was Metal Surfaces, Inc., Specification 



Fig. A-4. Preparation of plate for receiving thermocouple plug 

SP 10 1.) 1 2 After this step, all parts were handled 
with clean white gloves. 

(4) The parts were assembled for brazing. The thermo- 
couple plug was held down firmly against the bot- 
tom of the hole in the plate and locked in position 
by making a series of punch indentations in the 
plate around the edge of the hole. A finishing tap 
on the top of the plug with a light hammer, using 
a fixture to protect the wires, assured that the plug 
was firmly seated. 

The cleaned copper feed-tube was inserted in its 
socket (see Fig. A-4), and the brazing wire or slurry 

1B Metal Surfaces, Inc., Bell Gardens, Calif. 
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was placed in the tube. 13 (When slurry is used, it 
should be needle-injected so as to fill from the 
bottom to the top of the feed-hole and tube, 
thus eliminating the possibility of bubbles in this 
passage.) 

(5) The assembled parts were placed in a hydrogen- 
atmosphere furnace for brazing. During this oper- 
ation, the 0.127-cm (0.005-in, ) thermocouple lead 
wires extending from the back of the thermocouple 
plug were protected from damage, 

(6) After the brazing operation was completed, the 
woven fiber glass sleeve insulation was placed over 
the thermocouple lead wires and the final incre- 
ment of insulating ceramic cement was placed in 
the top of the thermocouple plug, as previously 
described. (At this time the protruding brazing- 
alloy feed-tube can be snipped off.) When the ther- 
mocouple lead wires have been connected to the 
appropriate terminals and/or cold junction and 
instrumentation lead wires and covered to protect 
them from damage, the thermocouple plug instal- 
lation can be considered complete. 

A few comments on the final design of the thermocouple 
plug should be added here, lest the reader be tempted 
to take shortcuts and be forced to retrace ground that 
was covered in the present investigation. 

The laser-fused closed tip is a most important feature 
of the thermocouple plug developed here. Numerous 
earlier attempts to incorporate prewelded thermocouple 
junctions into a plug configuration failed because out- 
gassing from ceramic insulating material destroyed the 
integrity of the brazed joints, or were impractical because 
of the several temperature levels of brazing operations 
which were required. The laser-fuzing of the thermo- 
couple junction and the tip of the thermocouple plug 
provides a high-temperature closure that is made under 
conditions where any possible outgassing from the single 
piece of double-bore ceramic tubing is harmlessly vented 
away from the weld area. The final brazing of the thermo- 


iS Two different brazing alloys, having different melting tempera- 
tures, were characterized during the brazing procedure develop- 
ment tests. The choice of alloy depends on the anticipated 
temperature to which the brazed joint will be exposed during 
test: BT alloy (72% Ag, 28% Cu) is supplied in wire form and 
flows at 1052 K (1435°F). Optimum diametral clearance for this 
alloy appears to be 0.00127 to 0.00254 cm (0.0005 to 0.0010 in.). 
Premabraze 130 alloy (82% Au, 18% Ni) is used as a slurry 
and flows at 1222 K ( 1740° F) . Optimum diametral clearance 
for this alloy appears to be 0.00254 to 0.00508 cm (0.001 to 
0.002 in.). 


couple plug into the recess in the plate is then not limited 
by the melting point of any previous brazing operation; 
the brazing alloy can be chosen freely to obtain the 
desired melting temperature and bonding characteristics. 
Further, any outgassing from the ceramic cement within 
the plug at the time of the brazing operation is sealed 
off from the braze and must vent to the back through 
the bore in the center of the plug. 

During the evolution of the brazing procedures, a 
number of types of fixtures for holding the plug firmly 
against the bottom of the hole in the plate were con- 
sidered and tried. None of the external fixtures was satis- 
factory, both because of bulk and mechanical complexity, 
and because loading derived from elastic deflection could 
not be depended upon at brazing temperatures. The 
assembly technique finally adopted was to mechanically 
wedge the plug in place by center-punching around the 
top edge of the hole in the plate; this was simple and 
reliable. 

Several methods of supplying brazing alloy to the 
joint area were tried before the external feed-tube was 
adopted: no other method provided the versatility of 
using wire or slurry brazing alloy, provided complete 
coverage of the plug-plate interface surfaces, and kept 
the reservoir cavity and thermal disturbance region as 
far from the thermocouple junction area. Test-specimen 
braze joints obtained using the procedures described 
here are shown in Figs. A-5-A-7. In these figures, the 
light-colored band at the interface between the plug and 
the body indicates a change in metallographic structure 
where the brazing alloy has diffused into the copper 
parts; the band width is many times the initial clearance 
between the parts. 

Satisfactory brazing assembly joints were demonstrated 
with two brazing alloys and over a range of diametral 
clearances (brazing alloy joint thickness). As a practical 
matter, it is desirable to strive for the minimum workable 
joint thickness, if only to assure that possible cocking of 
the plug in the hole is minimized. Some control over 
cocking can also he exercised during the punch-wedging 
operation by punching lightly first on one side and then 
on the other side of the plug. 

Absolute cleanliness of the brazing surfaces and alloy 
is essential if a complete bond is to be obtained. A clean- 
ing process at least as good as the one referenced must 
be used. Copper must be used for the brazing alloy feed 
tube; it was found that stainless steel contaminated the 
braze alloy and caused an imperfect bond. 
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Fig. A-5. Bond obtained with 0.635-cm (0.250-in.) diameter tellurium-copper plug brazed into hole in OFHC copper block; 
0.0025-cm (0.001-in.) diametral clearance, BT braze alloy (72% silver-28% copper) 
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Fig. A-6. Bond obtained with 0.635-cm (0.250-in.) diameter tellurium-copper plug brazed into hole in OFHC copper block; 
0.0051-cm (0.002-in.) diametral clearance, BT braze alloy (72% silver-28% copper) 
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Fig. A-7. Bond obtained with 0.635-cm (0.250-in.) diameter tellurium-copper plug brazed into hole in OFHC copper block; 
0.0051-cm (0.002-in.) diametral clearance, Premabraze 130 alloy (82% gold-18% nickel) 
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Appendix B 

Thermal Properties of Useful Metals 


A fundamental assumption in the transient heat transfer 
analysis and in the resultant curves presented in the body 
of this report is that the thermal properties of the plate 
material are constant. A further constraint on the plate 
materials which can be used in practice is that the level 
of the thermal conductivity and the thermal diffusivity 
must be high; this ensures that the thermal response is 
felt throughout the plate before the temperature of the 
heated surface exceeds its limiting value. 


Nickel is another material from which heat transfer 
apparatus can be fabricated. While the melting point of 
nickel is considerably higher than that of copper {1717 K, 
or 2630°F, vs 1358 K, or 1981 °F), the strength of nickel 
drops rapidly at temperatures above about 616 K (650 °F), 
and at 1200 K (1700 °F) the tensile strength is only 10% 
of the strength at room temperature. Nickel has a mag- 
netic transformation point at about 880 K (1125°F), and 
in this region the thermal properties are perturbed. 


The number of materials of construction which can 
be used to fabricate transient heat transfer devices and 
which will withstand fairly high temperatures, have high 
thermal conductivity, and have nearly constant thermal 
properties over a wide temperature range is relatively 
small. The thermal properties of three materials which 
have been used for rocket motor combustion chambers 
and nozzles are tabulated here for reference. 


Copper appears to be the most satisfactory material 
for the applications discussed herein. The thermal con- 
ductivity is higher than that of any other commonly 
available material, and the thermal diffusivity decreases 
by only 20% over a temperature range from ambient to 
about 1258 K (1800°F). The principal limitation of cop- 
per is in its low strength and its low “softening” tempera- 
ture. The yield strength of OFHC copper starts to decline 
rapidly when the temperature exceeds about 811 K 
(1000°F). For this reason, the surface temperature of a 
copper heat transfer apparatus should not be allowed to 
exceed a temperature of 811 K, except for a very short 
time. 


Mild steel is a readily available material, but in use 
has two limitations: the thermal conductivity and diffu- 
sivity are an order of magnitude lower than that of 
copper, and there is a phase change which causes a dis- 
continuity in thermal properties at a temperature of 
998 K (1337°F), which effectively limits the maximum 
temperature of operation. Over the operating temperature 
range, the thermal diffusivity changes by a factor of more 
than two. For these reasons, mild steel experimental 
apparatus can only be used when the heat flux level is 
low (say, up to 196 W/cm 2 , or 1.20 Btu/in. 2 s) and there 
will always be more uncertainty in the results than when 
the experimental apparatus is made of copper. 


The thermal properties of copper, mild steel, and nickel 
as collected from the literature are given in Table B-l. 
The data in this table show that the variation in thermal 
properties is much less for copper than for mild steel 
or nickel. 


A comparison of the transient heat flux capabilities of 
copper, mild steel, and nickel is given in Table B-2. To 
achieve a common basis for the comparison, it is assumed 
that the heated surface of each material reaches its maxi- 
mum allowable working temperature in 8 s, a time which 
should be large compared to any initial buildup period 
of the hot gas temperature and heat transfer coefficient 
in an experimental apparatus, and a time which is well 
within the response capability of ordinary data-recording 
equipment. It is further assumed that the 8-s time cor- 
responds to Fo = 1.00, thus assuring that the data are 
obtained into a region where the methods of analysis 
described here are applicable. The comparison is made 
for a flat plate having an initial temperature of 
T 0 — 256 K (0°F), and the effect of a wide range of 
hot-gas temperatures is considered. The first step in the 
comparison is to determine the characteristic plate thick- 
ness for each material. The definition of the Fourier Num- 
ber can be rearranged to give the plate thickness as 



The temperature ratio at the heated surface of the 
plate at Fo = 1.00 is evaluated as 



using the allowable maximum temperature for the mate- 
rial and the assumed gas temperature, and the Biot 
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Table B-l. Thermal properties of useful metals 




T 



p 



c o 


k 

a 

k 

Material 














K 

° F 

g/cm s 

lbm/in. 3 

joule/g K 

Btu/lbm °F 

W/cm 2 K 

Btu/in. s °F 

cm 2 /s 

in. 2 /s 


478 

400 

7.86 

0.284 

0.529 

0.126 

0.192 

6.54 X 10-« 

0.118 

.01826 


589 

600 





0.575 

0.137 

0.179 

6.08 X 10-4 

0.101 

.01562 


700 

800 





0.630 

0.150 

0.164 

5.59 X 10-4 

0.0847 

.01312 


811 

1000 





0.714 

0.170 

0.149 

5.08 X 10-4 

0.0679 

.01052 


922 

1200 

\ 


1 


0.823 

0.196 

0.132 

4.50 X 10-4 

0.0523 

.0081 


389 

240 

8.93 

0.3228 

0.397 

.0945 

1.48 

.005042 

1.068 

,1653 


500 

440 

8.88 

.3212 

0.408 

,0971 

1.46 

.004975 

1.029 

.1595 


611 

640 

8.83 

.3193 

0.419 

.0998 

1.45 

.004925 

0.997 

.1545 


722 

840 

8.77 

.3170 

0.430 

.1024 

1.44 

.004883 

0.970 

.1504 


834 

1040 

8.70 

'.3145 

0.441 

.1050 

1.43 

.004858 

0.949 

.1471 


945 

1240 

8.62 

.3117 

0.452 

.1077 

1.42 

.004842 

0.931 

.1443 


1057 

1440 

8.55 

.3089 

0.464 

.1103 

1.42 

.004817 

0.913 

.1415 


1168 

1640 

8.48 

.3061 

0.474 

.1129 

1.41 

.004808 

0.895 

.1389 


1279 

1840 

8.39 

.3033 

0.485 

.1157 

1.41 

.004796 

0.880 

.1365 


111 

-260 

8.88 

0.321 

0.252 

.0600 

0.453 

.001551 

0.520 

0.0807 


139 

-210 





0.311 

.0740 

0.409 

.001390 

0.378 

0.0586 


208 

-85 





0.395 

.0940 

0.320 

.001088 

0.233 

0.0361 


278 

+40 





0.433 

0.1030 

0.269 

.000915 

0.179 

0.0277 


417 

290 





0.494 

0.1175 

0.219 

7.46 X 10-‘ 

0.127 

0.01975 


556 

540 





0.571 

0.136 

0.188 

6.41 

0.095 

.01469 


625 

665 





0.609 

0.145 

0.179 

6.10 

0.085 

.01310 


695 

790 





0.534 

0.127 

0.185 

6.30 

0.100 

.01548 


834 

1040 





0.525 

0.1250 

0.202 

6.86 

0.110 

.01710 


972 

1290 





0.546 

0.1300 

0.218 

7.41 

0.115 

.01776 


1111 

1540 





0.566 

0.1350 

0.231 

7.87 

0.117 

.01820 


1250 

1790 





0.596 

0.1420 

0.247 

8.41 

0.119 

.01846 


1390 

2040 





0.650 

0.1500 

0.264 

8.96 

0.122 

0.01882 


1528 

2290 





0.664 

0.1580 

0.278 

9.45 

0.115 

.01865 


1668 

2540 





0.697 

0.1660 

0.293 

9.96 

0.121 

.01870 


2222 

3540 





0.840 

0.200 

0.357 

12.15 

0.122 

.01894 


Number corresponding to this temperature ratio at 
Fo ~ 1.00 is determined by interpolation on Fig, 9. From 
the Biot Number and from the thickness and thermal 
conductivity of the plate, the gas heat transfer coefficient 
is determined, and then the surface heat transfer rate is 
calculated. 

The comparison shows that the limiting heat trans- 
fer rate for copper (approximately 490 W/cm 2 , or 3.0 
Btu/in. 2 s) is twice that which can be absorbed by mild 
steel or nickel without overheating the inner surface of 
the material. 


For all three metals the maximum surface temperature 
is limited to 811 to 922 K (1000 to 1200°F), either by 
structural weakening or by a discontinuity in thermal 
properties which invalidates one of the assumptions of 
the analysis. 

Since copper can be used to nearly as high a wall 
temperature as the other metals, has a high heat flux 
capability, and has much less variation of thermal prop- 
erties, it is clear that copper is the preferred metal for 
fabrication of transient heat transfer measurement experi- 
mental apparatus. 
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Table B-2. Comparison of transient heat flux apparatus materials 11 



L 

T i 


1 1 
b b 



h 

( T i- 

T ) 

Q 


Material 

cm in. 

K 

°F 

Bi 

W/cin 2 K 

Btu/ 
in. 2 s °F 

K 

°F 

W/cm 2 

Btu/ 
in. 2 s 

Mild steel ( 1020 alloy ) 

.820 0.323 

1368 

2000 

.60 

.86 

.439 

.001492 

445 

800 

195 

1.195 

T w (max) = 922 K ( 1200°F) 
k = 0.165 W/cm* K (0.00056 Btu/in, s °F) 
a = 0.0838 cm 2 /s (0.013 in.Vs) 


3035 

5000 

.24 

.22 

.113 

.000382 

211 

3800 

237 

1.45 

Copper 

2.83 1.1135 

1368 

2000 

.50 

.61 

.790 

.002685 

556 

1000 

437 

2.68 

T w (max) = 811 K ( 1000°F) 


1923 

3000 

.33 

.33 

.427 

.001452 

1112 

2000 

475 

2.91 

k = 1.44 W/cm 2 K (0.0049 Btu/in. s °F) 


2480 

4000 

.25 

.23 

.298 

.001012 

1667 

3000 

496 

3.04 

a = 1.00 cm 2 /s (0.155 in.Vs) 


3035 

5000 

.20 

.18 

.233 

.000792 

2224 

4000 

516 

3.16 

Nickel 

.940 0.370 

1368 

2000 

.60 

.86 

.505 

.001720 

445 

800 

225 

1.377 

T w ( max ) = 922 K ( 1200° F ) 
k = 0.218 W/cm 2 K (0.00074 Btu/in. s °F) 
a = 0.110 cm 2 /s (0.0171 in.Vs) 


3035 

5000 

.24 

.22 

.129 

.000440 

211 

3800 

273 

1.672 

Comparison made on the basis of taking 8 s to attain Fo — 1.00 and 
T w — T max for the material, assuming T 0 = 0°F (flat plate). 



Equations used ; 

E*« 

II 

i 1 

Vi 

h 

_ Bik 
~L~ 


q = h(T s 

- TJ 


ui 

w 



Appendix C 
Sample Problem 


The analysis and procedures described in this report 
are applied here to a set of data obtained at the Jet 
Propulsion Laboratory during a series of tests conducted 
to determine both performance and heat transfer for a 
given rocket propellant, in this case Flox (88% F 2 , 
12% Oo-monomethylliydrazine, in a rocket motor de- 
signed to produce about 267 N (600 lbf) thrust. 

The thrust chamber and nozzle configuration are shown 
in Fig. C-l, which also shows the location of the thermo- 
couples. The thrust chamber is made of OFHC copper, 
with a wall thickness of approximately 2.54 cm (1 in.) 
(as shown to be optimum in Table B-2). Figure C-2 and 
Table C-l show the actual dimensions associated with the 
thermocouple plug installation. Note that the thermo- 
couple junction location is assumed to be 0.00889 cm 
(0.0035 in.) above the bottom of the hole. In the conical 
nozzle section, the equivalent cylindrical-plate radius of 


curvature is taken as the distance to the axis normal to 
the surface of the nozzle (shown as R* in Fig. C-2 and 
Table C-l). 

Data obtained during one firing of the rocket engine 
are plotted on Figs. C-3 and C-4. Figure C-3 shows the 
combustion chamber pressure as a function of time; it 
can be seen that the rise in pressure at the start of the 
test is quite rapid, which can be taken as an indication 
that the analytical requirement for sudden application 
of hot gas and convective heat flux is adequately satisfied. 
On the time scale on which the data were taken, the 
effective starting time for fully developed convective heat 
transfer appears to be approximately t 0 = 0.55 s. Fig- 
ure C-4 shows the temperature response vs time of a 
thermocouple plug located at the throat plane (Station ® 
on Fig. C-2); the initial temperature is seen to be 
T 0 = 302° K (84°F), and the midwall temperature at the 
end of the test is 545 °K (520° F). 


RING 

CLAMP 


OXIDIZER IN 


FUEL IN 


BOUNDARY LAYER 
COOLANT IN 


INJECTOR FACE 
COOLANT IN 


KISTLER PRESSURE 
TRANSDUCER 


ACOUSTIC CAVITY 
FILLER BLOCK 


THERMOCOUPLE 

PLUGS 



NOZZLE SECTION 


Fig. C-l. Thrust chamber configuration 
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Table C-l. Thermocouple plug installation dimensions (see Fig. C-2) 


Station 


L 




h 

X 


L/Ri 

x/L 

cm 

in. 

cm 

in. 

cm 

in. 

cm 

in. 

1 

2.73 

1.078 

3.43 

1.350 

1.53 


1.525 

0.600 

.798 

.557 

2 

2.81 

1.108 

2.97 

1.170 

1.57 

0.619 

1.565 

0.616 

.946 

.557 

3 

2.77 

1.091 

2.54 


1.57 

0.619 

1.565 

0.616 

1.091 

.565 



Table C-2. Mean thermal properties of copper used in 
sample problem data analysis (T = 340° F) 


Property 

Value 

P 

8.91 g/cm 3 (0.3220 lbm/in. 3 ) 

c p 

0.418 joule/g K (0.0958 Btu/lbm °F) 

k 

1.472 W/cm- K (0.00501 Btu/in. s °F) 

k 

a ~ 

1.408 cmVs (0.1624 in. 2 s) 


temperature rise at the test times corresponding to the 
selected Fourier Numbers. These quantities are also 
entered in Table C-3. This completes the collection and 
placing into usable form of the experimental data. The 
remainder of Table C-3 illustrates how these experimental 
data are used, in conjunction with the prepared ana- 
lytically based charts, to determine the effective hot gas 
temperature and heat transfer coefficient. 


Fig. C-2. Thermocouple plug installation in nozzle 

Following the steps outlined in Procedure II, the time- 
space mean wall temperature is estimated as 444°K 
(340 °F). The thermal properties of copper at this temper- 
ature are given in Table C-2. 

The time on the test-time scale corresponding to a given 
Fourier Number for convective heat transfer can be 
determined from the relationship 

(t -t 0 ) = — Fo 
a 

The value of L given for the throat-plane thermocouple 
plug in Table C-l, and the value of thermal diffusivity 
a given in Table C-2, are used to calculate the test time 
corresponding to Fo = 0.60, 0.80, and 1.00, as shown in 
lines 1 and 2 of Table C-3. The thermocouple plug tem- 
perature response record. Fig. C-4, is then used to deter- 
mine the corresponding temperature rise and the rate of 


The heat flux is calculated using the heat flux ratio 
from Fig. 30, as shown in Step 3 of Procedure II and in 
lines 6 and 7 of Table C-3. The heat flux determined 
in this manner is essentially independent of Fo and Bi 
and depends mostly on the measurement of the rate of 
temperature rise at the x/L = 0.55 position within the 
wall. 

The next step in the analysis of the test data is to 
calculate the quantity 

~2a (T x — T n ) ~ 

U / dT x \ 


at each Fourier Number, and to compare this value with 
the analytically calculated value at the same value of Fo, 
so that the matching value of Bi can be determined. This 
is done for the data at each of the three values of Fo in 
lines 8 and 9 of Table C-3. 
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TEMPERATURE, °F CHAMBER PRESSURE, psia 



E 


TIME, s 


Fig. C-3. Combustion chamber pressure, rocket engine test 
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TEMPERATURE, K CHAMBER PRESSURE 


Table C-3. Calculation of heat flux, hot-gas temperature and gas heat transfer coefficient 
from thermocouple plug test data at selected Fourier Numbers 3 


Line 

No. 


Parameter 


1 ~*o) = ^r F °’ S 

2 *tesl* s 

3 =0.565, test’ K { F) 

4 (T.-T 0 ) te , t .K(“F) 


\ dt /test 


P C p L ( 


/ dT. 
\~dt 


) 

/ J chart 

r r ( dT *W 

g = L pCp v^/J 


K/s ( °F/s) 

1 

from Fig. 30) 


X 


test 

= o.sr.n 


P c v L 


©J 


, W/cm 2 


m 


chart 

I.X/L,) =0.55 


~2« (t x - r 0 ) ~ 
L 2 


\ dt J _ test 
9 Bi ( from figure number noted ) 


10 = —Bi, W/cm 2 K 

Lt 


(-*L\ 

Vin. s °F/ 


11 


12 


( from figure number noted ) 


{x/D-oxb 


[ (T*-T 0 ) l 

L<r f - r 0 )J 

| ( from figure number noted ) 

\_(T g — r o ) J (*/£,) = i.ou 


13 (T,-T 0 ),K(’F) 

14 T g , K ( °F) 

15 (r,-T 0 ),K('F) 

16 T u „ K ( °F) 

17 (T ff -TJ,K(°F) 

18 q = h(T p -TJ, W/cm^-j 


(S) 


Fo = 0.60 

Fo = 0.80 

Fo = 1.00 

4.405 

5.870 

7.340 

4.955 

6.420 

7.890 

456 (361.5) 

499 (438.5) 

541 (514.0) 

410 (277.5) 

452 (354.5) 

495 (430.0) 

29.7 (53.40) 

28.7 (51.62) 

28.3 (50.95) 

1.55 

1.55 

1.55 

470 (2.88) 

439 (2.69) 

435 (2.66) 

1.417 

1.870 

2.30 

0.80 (41) 

0.60 (42) 

0.43 (43) 

1.080 (0.00367) 

0.814 (0.00277) 

0.580 (0.001975) 

.220 (31) 

.230 (32) 

.220 (33) 

.379 ( 36) 

.354 (37) 

.304 (38) 

956 (1260) 

nil (1540) 

1340 (1955) 

1002 (1344) 

1160 (1624) 

1389 (2039) 

521 (477.5) 

558 ( 545.0) 

586 (595.0) 

569 (562) 

605 (629) 

632 ( 679) 

691 (783) 

809 (995) 

1011 (1360) 

469 (2.87) 

451 (2.76) 

438 (2.68) 


a Given data for Table C-3 calculations: 
L = 2.77 cm ( 1.091 in. ) 

L/Rj = 1-091 
x/L = 0.565 
t 0 = 0.55 s 
r 0 = 302 K (84°F) 


JPL TECHNICAL REPORT 32-1588 


57 


The gas heat transfer coefficient is now calculated from 
the Biot Number, using the equation 

7 k 
h =L Bi 

The values of gas heat transfer coefficient are shown on 
line 10 of Table C-3. 


The surface heat flux corresponding to the heat transfer 
coefficient and to the gas and wall temperatures deter- 
mined above is calculated from the basic equation 

9 ~ g — T w ) 

which can be written in terms of the measured and 
charted quantities as 


The next step is to use the prepared charts to look up 
the temperature ratio at x/L = 0.55 and at x/L = 1.00 
(the hot wall) for each value of Fo and its corresponding 
value of Bi; these quantities are shown on lines 11 and 
12 of Table C-3, 



Now the actual calculation of the gas temperature and 
the hot-wall temperature can proceed. The gas tempera- 
ture is calculated by dividing the measured temperature 
rise at the thermocouple plug thermocouple junction by 
the temperature ratio for the x/L = 0.55 position as 
determined from the chart. Thus 

(T* ~ Tp) test , (x/L) =0.55 

(T, - To) 

(T, - T 0 ) 

and 

To = [(T, - T 0 ) + T 0 ] 


This heat flux is given on line 18 of Table C-3, and the 
values are seen to be acceptably close to the primary 
values of the heat flux given on line 7 of this table. 

The results of this sample problem data analysis can 
be reviewed to evaluate the validity of the procedures: 
It is believed that the rocket motor operation was such 
as to establish steady state combustion almost instanta- 
neously at the start of the test, with constant gas tempera- 
ture and a heat transfer coefficient which decreases 
gradually as the wall temperature increases (see below). 
The properties of OFHC copper are relatively constant 
over the range of plate temperature covered by the test. 
Thus the assumptions basic to the analysis are fairly well 
satisfied, 


(To ~ T 0 ) = 


Jchar;, (x/L) .;0.S5 


Similarly, the hot-wall temperature is calculated by 
multiplying the measured temperature rise by the ratio 
of the temperature ratios at the hot-wall and the midwall 
positions. Thus 


(T» To) — (T x T 0 ) lest 

<*/£,) -O.Sfii 


'(Tx ~ T 0 ) 1 


L(T*-T 0 )J 

(x/ L) - l.oo 

~(Tx — T 0 y 


L (To - To) 

(X/L) =0.5 3 / 


and 


Tw = [(T w - To) 4- T 0 ] 

From these results, the driving temperature difference 
for heat transfer, (T g ~T w ), can be determined. The 
temperature calculations for the sample case are given 
in lines 13 through 17 of Table C-3. 


The heat flux calculations agree, as they should, and 
seem to be of the right magnitude for the propellant and 
motor configuration being tested. There is no reason to 
consider the heat flux data to be in error by more than 
the percentage error in the measurement of the rate of 
temperature rise and in the evaluation of the average 
physical properties of the material, 


At increasing values of the test time and Fourier Num- 
ber, the apparent Biot Number decreases rather than 
remaining constant. While this may be partly due to the 
real effect of the hot-gas heat transfer coefficient decreas- 
ing slightly as the wall temperature increases, due to 
changes in gas density and viscosity with change in the 
average film temperature, it is more likely associated 
with high uncertainty ratios for determination of Bi in 
the range of Fo and Bi corresponding to the test con- 
ditions: in the range of Fo = 1.0 and Bi — 0.50, the 
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percentage uncertainty in Bi as determined by the meth- 
odology described here is approximately five times the 
sum of the percentage uncertainties in the Fourier Num- 
ber and in the ratio of temperature rise to rate of 
temperature rise, Since the Bi uncertainty ratios decrease 
with increasing Fo, and since the accuracy of the evalu- 
ation of both Fo and of the ratio of temperature rise to 
rate of temperature rise should be expected to increase 
with increase in Fo, it is to be expected that the final 
Biot Number on line 7 of Table C-3 is a better value than 
the preceding values, although this final value may have 
an uncertainty of as much as dz 50%. 

The apparent Biot Number variation during the test, 
and its uncertainty, are transmitted directly to the gas 
heat transfer coefficient, and, through the intermediate 
temperature ratio charts, this effect is felt in the deter- 
mination of the gas temperature and the hot wall tem- 
perature. 


The apparent gas temperature is seen to increase during 
the test, though even the final calculated value of 1388 K 
(2039 °F) is lower than the steady combustion tempera- 
ture which could be expected from the propellant used 
for the test. 

The conclusion to be drawn from die analysis of the 
test data as given in Table C-3 is that the heat flux is 
probably determined with good accuracy, but that there 
is considerable uncertainty in the gas heat transfer co- 
efficient and the hot-gas temperature. The fact that the 
gas temperature seems to be approaching a reasonable 
anticipated value as the test duration increases indicates 
that an increased test duration would have been desirable, 
widiin the limit imposed by the maximum allowable 
increase in inner wall temperature. This is consistent 
with the analytical study, which showed that Biot Num- 
ber determination accuracy at Fo =? 2 is much improved 
over a similar determination at Fo = 1. 
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